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Abstract 

THE  OPTICAL  PROPERTIES  OF  YOUNG  SEA  ICE 
By  Donald  Kole  Perovich 


Eight  laboratory  experiments  were  performed  to  determine  the 
optical  properties  of  young  salt  ice  and  to  examine  correlations  between 
the  optical  properties  and  the  state  of  the  ice.  Ice  grown  at  different 
temperatures  (-10,  -20,  -30,  and  -37°C)  and  from  water  of  different 
salinities  (0,  16,  and  31  °/00 )  was  investigated.  The  experiments  were 
conducted  using  a  tank  system  which  was  designed  to  grow  ice  similar  to 
that  found  in  nature,  allow  an  accurate  determination  of  the  state  and 
structure  of  the  ice,  and  permit  in  situ  optical  measurements  to  be  made. 
Measurements  of  incident,  reflected,  and  transmitted  irradiances  were 
used  in  conjunction  with  a  modified  Dunkle  and  Bevans  photometric  model 
to  determine  spectral  albedos  and  extinction  coefficients.  The  thin 
cold  ice  of  these  experiments  had  albedos  which  were  comparable  to  the 
values  for  the  thicker  warmer  ice  examined  by  previous  researchers; 
however,  extinction  coefficients  were  1.5  to  15  times  greater.  Qualita¬ 
tive  relationships  between  the  optical  properties  and  the  physical  state 
of  the  ice  were  observed.  As  the  ice  temperature  decreased  (and  the 
brine  volume  increased)  both  albedo  and  extinction  coefficient  increased; 
when  the  ice  temperature  dropped  below  the  eutectic  point  the  increase 
was  drastic.  The  dependence  of  albedo  and  extinction  coefficient  on  the 
brine  content  of  the  ice  was  found  to  be  complex  with  both  the  brine 


volume  and  its  vertical  distribution  being  significant.  Variations  in 


ii'i>  salinity  ovi'i  the  range  A  °/00  to  1A  l'/tn,  «.( i d  not  influence  the 
opt  leal  proport (os.  Increases  (n  albedo  anil  oxt loot  ton  coefficient 
woro  primarily  a  result  ol  changes  In  the  ire  structure  which  enhanced 
seat t  or ing. 

A  four  stream  discrete  ordinates  photometric  model  was  developed 
to  treat  the  case  ot  a  floating  ice  slab.  The  model  Included  both 
anisotropic  scattering  and  refraction  at  the  boundaries.  The  effects 
on  albedo  and  transmittance  of  variations  in  such  model  parameters  as 
the  single  scattering  albedo  and  the  phase  function  were  investigated. 
Using  one  and  two  layer  models,  theoretical  albedos  and  t ransmt t t ances 
were  compared  to  experimental  values.  The  four  stream  model  was  further 
applied  to  investigate  such  topics  as  the  depth  dependence  of  irradlance, 
cases  of  direct  beam  incident  radiation,  the  spectral  dependence  of 
albedo  and  transmittance,  and  the  effects  of  surface  layers.  It  was 
suggested  that  future  laboratory  research  be  oriented  towards  determina¬ 
tion  of  single  scattering  albedos  and  phase  functions  for  an  extensive 
range  of  ice  types. 
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CHAPTER  ONE 


BACKGROUND 


1.1  INTRODUCTION 

The  overall  energy  balance  of  the  polar  regions  is  greatly 
influenced  by  the  reflection,  absorption,  and  transmissions  of  short¬ 
wave  radiation  by  sea  ice.  The  amount  and  spectral  distribution  of 
light  penetrating  into  and  under  the  ice  is  biologically  important  in 
determining  primary  productivity.  Thermodynamically  radiation  penetra¬ 
tion  is  of  interest  since  heat  absorbed  in  the  ice  or  upper  ocean  can 
retard  freezing  or  contribute  to  melting,  thus  affecting  the  equilibrium 
thickness  of  the  ice. 

The  ice  cover  of  the  Arctic  Ocean  is  both  vertically  and 
horizontally  complex.  Vertically,  the  ice  exhibits  variations  in 
crystal  structure,  brine  volume,  temperature,  and  bubble  density. 
Horizontally,  thickness  can  vary  from  open  water  to  thick  pressure  ice 
over  distances  of  only  a  few  meters;  surface  conditions  are  diverse 
including  melt  ponds,  bare  ice,  and  snow  covered  ice. 

Optically,  sea  ice  is  a  complex  medium.  Unlike  its  constituent 
parts  of  brine  and  fresh  ice,  sea  ice  is  dominated  optically  by 
scattering  rather  than  absorption.  The  intricate  and  highly  variable 
structure  of  sea  ice  consisting  of  air  bubbles,  ice  platelets,  and  brine 
pockets  and  channels  considerably  complicates  the  processes  of  radiative 
transfer  within  the  medium. 

The  majority  of  the  previous  research  on  the  optical  properties 
of  sea  ice  has  been  concerned  with  thick  ice.  Both  spectral  and 


wavelength  integrated  albedos  and  extinction  coefficients  have  been 
determined  for  first-year  and  multiyear  ice  under  a  variety  of  field 
conditions.  However,  due  to  the  physical  hazards  involved  in  field 
measurements,  the  optical  properties  of  young  sea  ice  are  virtually 
unknown.  Although  such  thin  ice  accounts  for  only  a  small  percentage 
of  the  ice  cover,  it  transmits  large  amounts  of  energy  to  the  ocean 
and  hence  its  influence  is  important  to  the  regional  heat  budget 
(Maykut,  1978). 

The  only  practical  way  to  study  the  optical  properties  of  young 
ice  is  through  laboratory  experiments  which  allow  environmental 
control  and  avoid  the  inherent  hazards  of  the  field.  This  report 
describes  the  results  of  a  laboratory  program  designed  to  investigate 
the  optical  properties  of  young  ice.  The  term  "optical  properties" 
refers  to  the  radiative  parameters  of  albedo  and  extinction  coefficient 
over  the  wavelength  interval  400  to  1000  nanometers.  In  addition  to 
the  determination  of  these  optical  properties,  their  dependence  on 
such  physical  parameters  as  growth  rate,  air  temperature,  salinity, 
and  brine  volume  was  examined.  An  attempt  was  also  made  to  relate, 
through  a  photographic  analysis,  variations  in  the  optical  properties 
to  changes  in  the  ice  structure. 

The  final  goal  of  this  research  was  to  acquire,  through  photometric 
modeling,  a  more  precise  description  of  and  greater  physical  insight 
into  the  radiative  processes  occurring  in  thin  sea  ice.  Being  unable 
to  treat  such  complications  as  anisotropic  scattering,  direct  beam 
incident  radiation,  thin  slabs,  and  refraction,  photometric  models 
previously  applied  to  sea  ice  are  inadequate  for  this  task.  A  more 
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powerful  model  based  on  the  discrete  ordinates  method  (Chandrasekhar, 
1950)  was  therefore  developed  for  this  purpose. 

1.2  PREVIOUS  EXPERIMENTAL  RESEARCH 
1.2.1  Albedos 

Due  to  the  importance  of  net  shortwave  radiation  in  heat  budget 
calculations,  measurements  of  albedos  have  long  been  a  routine  part  of 
data  collection  in  the  Arctic.  Typical  values  of  bulk  albedo  are 
.06-. 1  for  open  water,  .2-. 6  for  melt  ponds,  .5-. 7  for  melting  bare  ice, 
to  as  much  as  .9  for  freshly  fallen  snow.  Areally  averaged  albedos 
have  been  measured  from  high  towers  and  aircraft.  Hansen  (1961)  and 
Langleben  (1971)  have  formulated  empirical  relations  between  albedo, 
the  degree  of  puddling  and  the  amount  of  open  water.  Selected 
results  of  previous  measurements  of  wavelength  integrated  albedo  are 
summarized  in  Table  1.1. 

Grenfell  and  Maykut  (1978)  have  determined  spectral  albedos,  from 
400  to  1000  nm,  for  an  extensive  range  of  arctic  ice  and  snow  conditions. 
Their  results  indicate  a  wavelength  dependent  albedo  displaying  larger 
values  at  shorter  wavelengths.  Both  the  magnitude  and  the  shape  of  the 
albedo  curve  are  extremely  sensitive  to  the  presence  of  liquid  water  in 
the  surface  layer.  A  reduction  of  albedo,  much  more  pronounced  at  longer 
wavelengths,  occurs  when  surface  water  is  present.  Pond  covered  ice 
has  a  maximum  albedo  at  short  wavelengths  with  a  severe  decrease,  by 
more  than  a  factor  of  five,  from  500  to  800  nm,  eventually  tapering  at 
1000  nm  to  a  value  determined  only  by  specular  reflection.  For  bare 
ice  the  magnitude  of  the  albedo  increases  and  its  wavelengths  dependence 


sharply  declines. 
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1.2.2  Extinction  Coefficients 

While  considerably  fewer  in  number,  measurements  of  radiation 
extinction  in  thick  first-year  and  multiyear  ice  have  been  made. 
Untersteiner  (1961),  from  changes  in  ice  temperature  profiles,  calculated 
an  optical  bulk  extinction  coefficient  of  .017  cm- 1  for  multiyear  sea 
ice.  Using  photometric  cells  Thomas  (1963)  measured  hulk  extinction 
coefficients  of  .011  cm-1  for  pack  ice  and  .0219  cm-1  for  shore  fast 
ice.  Weller  and  Schwerdtfeger  (1967),  studying  131  cm  thick  antarctic 
sea  ice,  determined  values  of  bulk  extinction  coefficient  of  .012  cm-' 
and  .010  cm-'  using  data  from  a  thermopile  radiometer  and  temperature 
profiles  respectively. 

Spectral  measurements  of  light  penetration  through  first-year 
(Mnvkut  and  Grenfell,  1975)  and  multiyear  (Grenfell  and  Maykut ,  1978) 
sea  ice  were  taken  with  a  submersible  spectrophotometer .  The  wide  range 
of  ice  types  investigated  showed  minimum  extinction  coefficient  for  ice 
beneath  an  old  melt  pond;  increasing  for  first-year  blue  ice.  the 
interior  of  multiyear  white  ice  and  the  surface  granular  layer  of  multi¬ 
year  white  ice.  In  all  cases  the  spectral  extinction  coefficient  was 
fairly  constant  from  400  to  500  nm,  exhibiting  a  gradual  change  between 
500  and  600  nm,  then  a  dramatic  increase  beyond  600  nm. 

An  understanding  of  the  optical  properties  of  salt  water  and  fresh 
ice  is  critical  to  effectively  model  the  transfer  of  radiation  within 
sea  ice.  There  is  an  abundance  of  literature  concerning  the  extinction 
of  light  in  various  water  types  (Neumann  and  Pierson,  1966). 

Goodrich  (1970)  and  Hobbs  (1974)  provide  reviews  for  both  bubble  free 
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and  bubbly  fresh  Ice.  Spectral  extinction  coefficients,  from  400  to 
800  tun,  for  fresh  Ice  and  sea  Ice  are  summarized  In  Figure  1.1, 

1.2.3  Correlation  of  Optical  Properties  and  the  Physical  State  of  the  Ice 
To  date  there  has  been  little  research  done  in  attempting  to 
correlate  the  optical  properties  to  the  physical  state  of  sea  ice. 

Weller  (1972),  studying  young  sea  Ice,  observed  an  increase  of  bulk 
albedo  with  thickness.  Experiments  conducted  by  Davis  and  Munis  (1973), 
using  a  HeNe  laser  (633  nm)  and  ice  samples  frozen  in  a  55  gallon  drum. 
Indicated  that  extinction  coefficient  Increased  as  the  ice  salinity 
increased. 

Further  Investigations  of  the  relationship  between  salinity  and 
extinction  coefficient  were  conducted  by  Lane  (1975).  Several  small 
samples  of  various  salinities  were  grown  at  a  temperature  of  -20°C. 
Transmission  was  found  to  decrease  as  salinity  increased,  though  the 
dependence  weakened  for  salinity  greater  than  5  °/po.  However,  Lane's 
samples,  being  frozen  solid  from  the  outside  to  the  center,  are 
structurally  and  perhaps  optically  quite  different  from  natural  sea  ice. 


WAVELENGTH  (nm) 


.0001 


Figure  1.1:  Review  of  spectral  extinction  coefficients: 
(l)  bubble  free  fresh  ice  (Sauberer,  I'JSO);  (2)  first 
year  blue  ice;  (3)  interior  of  multiyear  white  ice; 
and  (A)  surface  granular  layer  of  multiyear  white  i 
( 2— A  from  Grenfell  and 


CHAPTER  TWO 


RADIATIVE  TRANSFER  MODELS 


2.1  INTRODUCTION 

Experiments  investigating  the  optical  properties  of  a  medium 
typically  measure  incident,  reflected,  and  transmitted  irradiances. 

From  these  data  the  albedo  and  percentage  of  radiation  transmitted  are 
readily  calculated,  but  additional  information  is  needed  for  a  more 
detailed  analysis.  To  calculate  and  understand  the  optical  properties 
and  processes  of  a  medium  from  measured  irradiances,  a  suitable 
photometric  model  describing  the  transfer  of  radiation  within  the 
medium  must  be  determined. 

Historically,  the  Bouguer-Lambert  law  has  been  the  simplest  and 
most  commonly  used  method  for  investigating  radiation  extinction  in 
snow,  fresh  ice,  sea  ice,  and  water.  While  its  simplicity  is  appealing, 
deficiencies  of  the  model  preclude  its  application  in  this  experiment. 

In  particular  the  Bouguer-Lambert  model  assumes  an  infinitely  thick 
medium,  which  is  invalid  when  examining  young  thin  sea  ice.  Two  photo¬ 
metric  models  are  used  to  analyze  the  experimental  results  of  this 
research,  the  Dunkle  and  Bevans  and  discrete  ordinates  methods. 

2.2  DUNKLE  AND  BEVANS  MODEL 

Dunkle  and  Bevans  (1956)  developed  a  photometric  model  to  analyze 
the  transmittance  and  reflectance  of  a  snow  cover.  This  technique  was 
modified  by  Grenfell  and  Maykut  (1978)  to  calculate  albedos  and 
extinction  coefficients  for  snow  and  sea  ice.  This  approach  was  also 
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selected  to  determine  extinction  coefficients  in  the  laboratory 
experiments . 

The  theory  assumes  an  optically  homogeneous  medium,  finite  in  depth, 
infinite  in  width  and  breadth,  and  isotropically  scattering;  the 
radiation  field  in  the  medium  is  assumed  to  be  diffuse  at  all  depths. 

In  the  present  case  the  model  is  adopted  to  treat  the  additional  compli¬ 
cations  of  specular  reflection  at  the  upper  surface  and  transmitted 
measurements  made  within,  as  well  as  beneath  the  ice.  The  geometry  of 
the  Dunkle  and  Bevans  model  for  this  experiment  is  illustrated  in 
Figure  2.1.  The  Dunkle  and  Bevans  model  assumes  two  streams,  an 
upwelling  and  a  downwelling  irradiance.  This  entails  the  solution  of 
two  simultaneous  differential  equations  (Dunkle  and  Bevans,  1956), 

dF^(z.A)  =  -k^Fv(z,A)dZ  -  r^F+(z,A)dz  +  r^F^(z,A)dz  (2.1) 

dFf(z,A)  =  kxF+(z,A)dz  -  rxF+(z,A)dz  +  rxFf(z,A)dz  (2.2) 

where  F^(z,A)  is  the  downwelling  spectral  irradiance,  F  (z,A)  is  the 
upwelling  spectral  irradiance,  k  is  the  spectral  absorption  coefficient, 

A 

rx  is  the  spectral  volume  reflectance  coefficient,  and  z  (increasing 
downward)  is  the  depth  within  the  ice.  The  boundary  conditions  for  this 
problem  are 

z  =  0  F+(0,A)  =  (1-R0)F0  -  R0F+(0,A)  (2.3) 

F+(H,A)  =  0 


z  =  H 


(2.4) 


Figure  2.1:  Geometry  of  Dunkle  and  Bevans  model. 
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where  Rp  is  the  specular  reflection  from  the  upper  surface,  F(1  is  the 
incident  downwelling  irradiance  at  a  given  wavelength,  and  H  is  the 
thickness  of  the  ice. 

For  a  homogeneous,  plane  parallel  slab  of  thickness  H,  with  the 
above  boundary  conditions  [Equations  (2.3)  and  (2.s)],  the  solutions  to 
Equations  (2.1)  and  (2.2)  are 

(l-R,)Fn  sinh[<  (H-z)+y] 

V*.X)  =  - - - * -  (2.5) 

sinh[ ]  -  R0sinh[<^H] 


Ff(z,A) 


where 


(1-R3)F0  sinh[<x(H-z) ] 
sinh[<xH+>]  -  R0sinhU.H] 


*  +  2k  r 

ka  ka  a 


(2.6) 


(2.7) 


is  the  extinction  coefficient  and  y  =  arcsinh  ( ic x / r x ) .  In  the  limiting 
case  of  very  thick  ice  (H  -♦  <*>)  Equations  (2.5)  and  (2.6)  reduce  to 


F+(z,A) 


1-R0 

1  -  R06x 


Fne 


-Kl  z 


Ff(z,A) 


(1-R0)6X 

1-R0PX 


_  -Kv  S 

F0e  X 


(2.8) 


(2.9) 


which  is  equivalent  to  the  Bouguer-Lambert  law  (where 

0A  =  t(rA  +  <X) l/2_<A]/rA) * 

Combining  Equations  (2.5)  and  (2.6)  relationships  are  written  in 
terms  of  the  measured  values  of  albedo  and  transmittance.  For  albedo, 
a(A) , 

R0F0+(1-R0)F  (0,A) 

«(A)  =  - r- - 

F0 


(2.10) 
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Evaluating  Equation  (2.6)  for  z  equal  to  zero. 


Ff(0,X)  = 


(l-R0)F0sinh(<AH) 
sinh^^H+y)  -  RgSinh^H) 


(2.11) 


and  substituting  into  Equation  (2.10), 


o(A)  =  Rn  + 


(l-Ro)2sinh(KxH) 
sinh(<xH+Y)  -  R0sinh(KAH) 


(2.12) 


The  transmittance,  T(z,X),  is  expressed  as 


F  (z,X)  ( 1 — Rf>  ) sinh [  k,  (H-z )+y ] 

T(z » X)  =  — - -  =  - - - ± -  .  (2.13) 

Fq  sinh(<xH+y)  -  R^sinh  k^H 

Equations  (2.12)  and  (2.13)  cannot  be  explicitly  solved  for  k^  and 
r^.  Instead  the  equations  are  numerically  inverted  using  a  two- 
dimensional  Newton-Raphson  scheme.  This  entails  the  iterative  solution 
of  the  simultaneous  equations 


(a‘f)i  ~  3^ 


Ak  - 


_3f 

3r 


Ar  =  0 


(2.14) 


(T-g). 


3k 


Ak 


_ 

3r 


Ar  =  0 


(2.15) 


where  a  and  f(K,r)  are  respectively  the  observed  and  calculated  albedo 
and  T  and  g(K,r)  are  the  observed  and  calculated  transmittance. 
Equations  (2.14)  and  (2.15)  are  expressed  in  matrix  form  by 


(2.16) 


n 


Initial  values  of  k  and  r  are  selected  using  the  Bouguer-Lambert  model 
and  substituted  into  Equation  (2.16).  Values  of  Ak  and  Ar  are  then 
determined,  <  and  r  adjusted,  and  the  process  repeated.  The  iteration 
continues  until  the  change  in  k  and  r  is  less  than  one  percent. 

Convergence  is  rapid,  typically  occurring  within  3  iterations. 

Using  the  Dunkle  and  Bevans  method  spectral  reflectances  and 
extinction  coefficients  are  readily  calculated  and  compared  to  results 
of  previous  researchers.  However  with  only  two  streams  the  IHinkle  and 
Bevans  model  provides  no  information  concerning  the  angular  distribution 
of  radiation  within  the  ice.  The  model  is  also  limited  by  the  assumptions 
of  diffuse  incident  radiation  and  isotropic  scattering.  There  is 
evidence  to  suggest  that  the  scattering  in  sea  ice  is  quite  anisotropic 
(Grenfell  and  Hedrick,  personal  communication).  These  constraints  limit 
the  applicability  of  the  Dunkle  and  Bevans  method.  In  order  to  accurately 
describe  the  radiative  transfer  processes  occurring  in  the  ice,  a  more 
general  photometric  model  is  needed. 

2.3  DISCRETE  ORDINATES  METHOD  (DOM) 

2.3.1  Background 

Chandrasekhar  (1950),  using  a  discrete  ordinates  method,  derived 
solutions  to  the  equation  of  radiative  transfer  in  a  plane  parallel, 
homogeneous ,  multiple  scattering  medium.  The  discrete  ordinates,  or 
streams,  refer  to  the  angles  at  which  radiances  are  determined.  Liou 
(1973,  1974)  developed  a  numerical  model,  based  on  Chandrasekhar's 
solutions,  to  investigate  radiative  transfer  within  cloudy  and  hazy 


atmospheres.  Analytic  solutions  were  det  ermined  tor  ttio  two  ami  tour 
stream  eases  (Lion, 

A  primary  interest  in  vertical  lrradianee  distributions  ami  the 
need  to  Investigate  a  large  number  of  cases  dictated  the  use  ot  a  four 
stream  model  In  this  research.  Numerical  experiments  hv  l.lou  (l‘*7Al 
established  the  effectiveness  of  the  four  stream  model  in  studies  ot 
lrradianee  distributions  in  turbid  media.  A  two  stream  model  was 
unacceptable  since  errors  became  prohibitively  large  for  optically 
thin  layers.  Higher  order  models,  while  offering  greater  angular 
resolution  of  radiance,  were  complicated  by  numerical  instabilities  and 
iterative  solution  schemes. 

The  model  used  here  Is  based  on  the  four  stream  solution  discussed 
by  l.lou  (1974).  However  the  new  application  of  the  model  to  sea  ice 
Introduces  several  complications,  primarily  associated  with  boundary 
conditions  and  refractive  effects.  These  difficulties  are 
discussed  In  Section  2.3.2,  while  the  remainder  of  this  sect  ton  is 
limited  to  a  general  description  of  the  discrete  ordinates  method.  The 
geometry  of  the  four  stream  model  is  illustrated  in  Figure  2.2. 

From  Chandrasekhar  1 1  SO')  the  equation  of  radiative  transfer  tor 
plane  parallel  atmospheres  is 

U  sU-tlolLiii  _  l(t  ,p,T>  (2.17) 

dT 

where  1  is  the  radiance,  $  is  the  azimuth  angle,  w  Is  equal  to  cos  0, 
and  0  Is  the  inclination  to  the  vertical.  The  sign  convention  Is  that 
-u's  represent  downward  angles.  The  optical  depth  i  can  be  expressed 


INCIDENT  RADIATION 


Figure  2,2:  Ceometry  of  four  stream  model. 


in  terms  of  the  scattering  coefficient  o,  and  the  absorption  coefficient 


k. 


t  “  -  (o+k)z 


(2.18) 


The  minus  sign  is  a  result  of  the  custom  that  t  increases  downward  and 
7.  decreases  downward.  The  source  function  S  describes  the  contribution 
of  multiple  scattering  and  for  a  plane  parallel  scattering  atmosphere 
with  a  direct  incident  beam  is  expressed  as 


J  .-2n 


•S(T,u,i)  -  IpCu.'i.u'.'f'iUt.p’.'f'idvi'di'' 


-1  0 


(2.19) 


+  £E0p(u,tftn0,'I'0)e 


-i/p 


For  this  analysis  axial  symmetry  is  assumed  and  Equation  (2.17)  reduces 
to 

1 

U  ci  r^-T-  l‘i  =  I  ( i ,  vi )  -  “  f  p  (p  .  W  1 )  I  (T  , p  1  )dp  ' 
dr  -i 

(2.20) 

If  /  \  /  Vl,. 

-  -E0p(p,vi0)e  o 

where  E0  is  the  solar  radiance  and  Hq  is  the  cosine  of  the  solar  zenith 
angle,  defined  to  be  positive. 

The  scattering  processes  within  the  medium  are  Jointly  defined  by 
the  phase  function,  p(h,vi')t  and  the  single  scattering  albedo,  u>()  .  The 
phase  function  describes  the  angular  dependence  of  scattering  and 
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typically  Is  normalized  so  that  the  Integral  over  angle  Is  equal  to  one. 
Different  types  of  scattering  from  Isotropic  to  strongly  anisotropic  may 
he  represented  hv  the  proper  choice  of  a  phase  function.  Expressed  in 
terms  of  Legendre  polynomials, 

N 

p(p.p')  -  y  uePe(p)Pt(p')  (2.21) 

<-0 

where  N  is  three  for  a  four  stream  model.  A  major  limitation  of  the 
four  stream  model  Is  that  a  four  point  expansion  cannot  accurately 
describe  complex  phase  functions. 

The  single  scattering  albedo  <Ii(1  describes  the  relative  importance 
of  scattering  and  absorbing  processes. 


u>0 


o 

o+k 


(2.22) 


<I>0  ranges  from  zero  for  a  purely  absorbing  medium  to  one  for  a  purely 
scattering  medium.  The  higher  order  <1^  's  used  in  Equation  (2.21)  can  be 
obtained  from  the  expression 

.1 

I  p(x)Pj(x)dx  ,  (2.23) 

-1 

which  is  derived  from  Equation  (2.21)  exploiting  the  orthogonality  of 
the  Legendre  polynomials.  For  simple  phase  functions,  such  as  isotropic 
or  Rayleigh  scattering,  the  integral  in  Equation  (2.23)  can  be  evaluated 
exactly,  while  for  more  complex  phase  functions  numerical  integration 
schemes  must  be  used.  The  ui^’s  for  1  greater  than  zero,  are  normalized 
by  multiplying  the  results  of  Equation  (2.23)  by  <o() . 


The  total  solution  Is  the  superposition  of  Equations  (2.26)  and  (2.27). 

Following  the  development  of  Chandrasekhar,  Liou  (1974)  determined 
the  analytic  solution  for  a  four  stream  model.  After  lengthy  algebraic 
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manipulation  an  expression  for  radiance  as  a  function  of  depth  and  angle 
was  determined; 


U>.n()  »  y  LJW)(nt)e'kjT  +  Z(n1>e",/vi°  .  (2.28) 

J— 2 
.1  f  0 

For  a  four  stream  model  the  functions  w j ( V* 4 )  and  Z(u{)  expand  to 

3 

yy  ■  TSTkT  Yyyyyy 

J  0  _  a 


(2.29) 


V»il 


2  /  2 

PoE0(h,-ho)(P2-Po) 


z(vli) - -- — - — 1 — — — —  y*,ct(-- 

(l~kjhp)  (1-k^pp)  ^  k  W< 


-Mr  (u  )  . 

0  *  1 


(2.30) 


The  kj's  denote  the  eigenvalues  of  the  homogeneous  differential  equations, 
while  the  {.^'s  are  proportionality  constants  which  arise  solving  the 
associated  homogeneous  equations  and  can  be  calculated  from  the  recursion 
formula 


2i+l-<a 

CC+1W  “  X(M-l)  “  C+T^t-l^  (2.31) 

where  £.0(X)  -  1  and  f._;(X)  -  0. 

The  L^'s  in  Kquat ion  (2.28)  are  constants  of  integration  determined 
from  the  boundary  conditions.  The  system  of  four  equations  [Equation 
(2.28)1  and  four  unknowns  (Lj's)  is  expressed  in  matrix  form  and  solved 
using  a  standard  inversion  routine.  The  determination  of  the  proper 
boundary  conditions  representing  the  downwelling  radiances  at  the  upper 
boundary  and  the  upwelllng  radiances  at  the  lower  boundary  is  discussed 


in  Section  2.3.2. 
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To  determine  the  upward,  downward,  and  not  irradianees  at  any  level 
it  is  only  necessary  to  express  the  integral  of  radiances  as  a  sum  of 
Gaussian  quadratures  and  evaluate.  The  net  irradiance  is  expressed  by 

1 

Fuet(T)  “  •%(*)  -  F4(t)  ”  2tt|  I(j,p)udp  -  nu0E0e  ,  (2.32) 

-1 

expanding  the  integral 

2 

Fnet(T)  "  2n  ajKT.u^Pj  -  nlpEpe"’ /u0 ,  (2.33) 

J— 2 

and  correspondingly 

Ff(i)  »  lutajPjHx.pj)  +  a_„u,I  (t  ,u ,)  ]  (2.3-4) 

F4(t)  -  -2w  tajUjI(T,-pj)  +  a,vi,I(T,-p^)  ]-  wp^E^e  '^"o  .  (2.3S) 

2.3.2  Determination  of  Boundary  Conditions 

With  a  general  analytic  four  stream  model  developed  by  Liou  (1^74), 
the  remaining  theoretical  difficulties  are  to  determine  the  proper 
boundary  conditions  and  derive  expressions  for  albedo  and  transmittance 
for  the  case  of  a  floating  ice  slab.  This  is  done  both  neglecting  and 
considering  refraction. 
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Case  of  No  Refraction 

If  refraction  at  the  air-ice  and  Ice-water  interfaces  is  ignored, 
the  boundary  conditions  are 


K-hj.O)  =  l(-y2,0)  =  C  (upper  boundary)  (2.36) 


I(lJl»TN)  =  I(p2,tn)  =  0  (lower  boundary)  (2.37) 


where  is  the  optical  thickness  of  the  ice  slab.  For  isotropic 
incident  radiation,  representing  cloudy  sky  conditions,  C  is  equal  to 
one  and  there  is  no  direct  beam  component.  To  represent  clear  skies,  C 
is  assumed  to  be  zero  and  the  incident  radiation  field  is  defined  by  a 
direct  beam  of  radiance  EQ  at  an  angle  arc  cos  (u0). 

The  albedo  and  transmittance  are  determined  using 


F.(0) 

- -  (2.38) 

F0 


T 


w 


(2.39) 


Equations  (2.34)  and  (2.35)  are  used  to  calculate  F^(0)  and  F4(tn), 
where  F„  is  equal  to  n  for  isotropic  incident  radiation  and  u  nE  for  an 

o  K  0  0 

incident  direct  beam. 

For  a  snow  surface  (Barkstrom,  1972)  it  is  possible  to  consider 
refraction  solely  as  a  scattering  phenomenon.  However,  for  sea  ice,  it 
is  usually  necessary  to  calculate  refractive  effects  at  the  interfaces 
explicitly. 


Considering  Kofvaot ion 

When  light  passes  from  one  medium  to  another  of  different  index  of 

refraction,  Fresnel  reflection  and  refraction  occur.  Figure  2.3 

illustrates  the  elementary  optical  situation  of  a  ray  of  light  incident 

on  a  boundary  between  media  with  different  indices  of  refraction.  The 

angles  0,  0  ,  and  $  are  related  by 
K 

e  -  0R  (2.40) 

and  from  Snell’s  law 

nj  sin  6  »  n,  sin  ,  (2.41) 


where  nj  and  ng  are  the  respective  indices  of  refraction. 

It  is  also  necessary  to  determine  the  relationship  of  the  incident 
(Ig),  reflected  (IR),  and  transmitted  (Ig)  radiances.  1R  is  calculated 
from  the  Fresnel  equations  (Born  and  Wolfe,  1965)  to  be 


iR<e>  -  i„<e>f 


sln~  (O-if)  tan-  (d->f) 
sin- (0+^)  tan*(0+4>) 


(2.42) 


The  transmitted  radiance  is  obtained  from  energy  conservation. 


ij(*>  "  (V9> 


cos  B 
cos 


(2.43) 


For  a  floating  ice  slab,  refraction  can  occur  at  ice-water  and  ice- 
air  interfaces,  situations  which  are  illustrated  in  Figure  2.4.  Since 
it  is  assumed  that  there  is  no  upwelling  radiation  from  the  water, 
refraction  at  the  water-ice  interface  need  not  be  considered.  Due  to 
the  small  difference  in  the  value  of  n  for  water  and  ice,  refractive 
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effects  on  downwelling  radiation  in  the  ice  reaching  the  ice-water 
Interface  are  quite  small.  For  example,  from  Equation  (2.42)  for  unit 
Iq *  Ro  *0001  at  u;  and  .0026  at  p j .  The  influence  of  refraction  is 
negligible  at  the  lower  interface,  leaving  the  downwelling  irradiance  at 
and  the  lower  boundary  conditions  [Equation  (2.37)1  unchanged. 

At  the  upper  interface,  with  a  much  larger  index  of  refraction 
mismatch,  refraction  is  a  significant  factor.  Often  in  photometric 
modeling  an  Isotropic  incident  radiation  field  is  assumed,  however  the 
transformation  of  this  isotropic  field  by  refraction  at  the  air-ice 
interface  is  usually  oversimplified.  The  isotropic  incident  radiation 
field  in  the  air  (before  refraction)  is  defined  by 


'air'".’ 


(2.4M 


where  ua  ■  cos  6.  Determined  by  either  geometric  considerations  (Jerov, 
1976)  or  by  variable  transformation  (Papoulis,  1965),  the  redistribution 
of  the  radiation  field  in  Equation  (2.44)  by  refraction  is 

'lce“V>  '  W“.>  «•«> 

where  pe  ■  cos  4>,  and  ranges  from  1  to  .646.  For  isotropic  incident 
radiation,  Figure  2.5  displays  the  incident  [from  Equation  (2.44)1, 
reflected  [from  Equation  (2.42)1  and  refracted  [from  Equation  (2.45)] 
radiation  fields.  Before  refraction  the  radiance  is  one  at  all  angles. 


while  after  refraction 


ANGLE  (deg) 


Figure  2.5:  The  angular  distribution  of  isotropic  incident 
radiation  undergoing  refraction  =  49.5°'):  (1)  isotropic 
incident  radiation;  (2)  incident  radiation  after  the  removal 
of  specular  reflection;  and  (3)  radiation  distribution  in 
the  ice  after  refraction 
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Iice(_ui ,0)  =  0  (2’46) 

Ilce(-V*2>0)  =  ll67  •  (2.47) 

Substituting  the  values  from  Equations  (2.46)  and  (2.47)  into 
Equation  (2.35),  the  downwelling  irradiance  is  calculated  to  be  -1.17tt. 
However,  it  can  be  shown  analytically  that  the  downwelling  irradiance  in 
the  air  is  equal  to  -it.  The  erroneous  result  of  -1.17tt  is  caused  by  the 
inaccuracies  of  what  is  essentially  a  one  point  quadrature.  To  avoid 
this  and  conserve  energy,  it  was  necessary  to  adjust  the  value  of 
I^e(-P2  ,Q) .  The  downwelling  irradiance  after  refraction, 

.646 

F0lce  =  -2tt|  Ilce(n)dp  ,  (2.48) 

was  more  precisely  evaluated  using  a  multi-interval  Simpson's  rule  as 
.937tt.  The  irradiance  determined  by  Gaussian  quadratures  reduces  to 

F0lce  =  -2ira2y2llce(-p2,0)  .  (2.49) 

Setting  Equation  (2.49)  equal  to  .937  yields 

Iice(-h2.°)  =  1-564  ,  (2.50) 

which  is  used  rather  than  Equation  (2.47). 

The  upper  boundary  conditions  are  further  complicated  by  refraction 
and  reflection  of  the  upwelling  radiation.  Since  arccos(pi)  is  greater 


.’8 

than  tin1  critical  angle  -  49.70"),  l(m,0)  i  totally  rot  loot  oil  an. I 
from  t ho  Fresnel  equation  .022  of  l(u.>,0)  is  reflected.  t'omhining  those 
results  with  Kquat ion  (2.40)  and  (2.  Ml)  yields  the  tippet  houmlarv 
coinlltlons  for  ditfuse  incident  radiation. 


1  (-ii  t  .0)  -  1  (n  t  ,0) 


(.'.mi 


1 (-P  ,  ,0)  -  .022  l  ui ,  ,0)  +  0 


t .  . 


whore  0  ■*  1.504. 

II  the  inetilont  radiation  is  a  direct  hoam  K,(u(),  the  rot  i  acted 

hoam  l-’  *  ( u  ’ )  is  determined  from  liquations  (2.41),  (2.42),  and  (2.4  U. 

The  boundary  conditions  are  the  same  as  for  diffuse  incident  |  1'quat  ions 

(2.51)  and  (2.52)1  except  that  for  a  direct  beam  0  is  equal  to  ,-oro. 

While  refraction  does  not  change  the  expression  tot  t t ansmi t t unco 

liquation  (2.  I‘>)  |,  it  does  complicate  the  calculat  ion  ol  albedo.  I'o 

determine  transmittance  and  albedo,  K  ,  K.(t„l,  and  I-'  (01  must  ho 

i1  t  N  ait* 

calculated.  K  and  K.(t,,)  remain  t  he  same  as  I  lie  nonrel  t  act  i  vo  case, 

I'  f  N 

leaving  transmittance  unchanged.  lo'wever,  is  quite  dttterent. 

It  is  important  to  remember  that  the  albedo  is  the  t  at  to  ol  the  incident 
dowttwelling  to  the  tipwelling  irradiance  in  the  ait  ,  before  the  incident 
light  is  refracted  and  after  the  upwolllng  is  ret t acted.  I  ho  model 
calculates  tipwelling  radiances  at  angles  pj  and  u<  in  the  ice.  In  the 
case  of  no  t  el t  act  ion 


Fairt<0)  '  Ktcef(0) 


2T,(ntv,l 1  Ico0  !  *0>  4  V1:'! ce(,,2*0')  * 


n 
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But  when  considering  refraction,  with  1 internally  reflected. 
Equation  (2.53)  no  longer  is  applicable. 

It  was  decided  that  the  albedo  could  most  accurately  be  calculated 
by  allowing  ljce(bi,0)  and  ^(pp,0)  to  define  an  upwelling  radiation 
field  of  the  form, 


w  ith 


and 


I.  ^(p)  -  Ati  +  B 
i  ce 


‘ice1".’0’  -  WV°P 


(P,  ~  P?) 


B  =  Ilce(lll’°)  "  A,'l 


(2.54) 

(2.55) 

(2.5b) 


The  form  of  Equation  (2.54)  is  the  highest  order  possible  with  only  two 
points.  As  a  simple  consistency  check  of  Equation  (2.54),  the  end  points 
1(1,0)  and  1(0,0)  are  evaluated  to  insure  that  they  are  non-negative. 

The  upwelling  irradiance  emergent  from  the  ice  is  expressed  as 

F,irf(0)  "  2"/  Iice<M)(1“RF)udl'  •  (2.57) 

l'c 

The  cosine  of  the  critical  angle,  p  ,  is  the  lower  limit  of  integration 
since  light  from  zero  to  p^  is  internally  reflected.  Rp  is  a  correction 
for  the  Fresnel  reflection  of  the  emergent  radiation  and  is  assumed  to 
equal  .065,  which  is  the  Integrated  value  for  Fresnel  reflection  under 
Isotropic  conditions.  Equation  (2.57)  now  becomes 

F  .(0)  -  2n ( .935) I  (Ap  +  R)pdp 
air*  ' 


(2.58) 


)0 

Integrating  gives 

F  .  .(0)  =  5.87  (.243  A  +  .291  8)  .  (2.59) 

air  t 

The  upwelling  irradiance  in  the  air  consists  of  F^rf(0)  plus  a 
contribution  from  Fresnel  reflection  of  the  incident  radiation. 


Fttotal  (0) 


F  ,(0) 
air  t 


+  R  ( 


(2.80) 


where  R0  «  .065n  for  isotropic  incident  radiation  and  Rp  =  u ^  •>  1  ^ ( u  ( ) 
for  an  Incident  direct  beam,  with  lu  determined  hv  Kquation  (2.42).  Bv 
substituting  for  F^(0)  in  Equation  (2.38)  the  albedo,  consider¬ 

ing  refraction,  is  calculated. 

Based  on  this  solution  to  the  equation  of  radiative  transfer  a 
system  of  computer  programs  was  written  to  investigate  the  optical 
processes  of  a  floating  ice  slab.  Model  input  parameters  are  optical 
thickness,  phase  function,  and  single  scattering  albedo.  The  four  stream 
model  calculated  albedo,  transmittance,  and  radiances  and  irradiances  at 
any  depth  within  the  ice.  Direct  or  diffuse  incident  radiation  is 
introduced  through  the  boundary  conditions.  The  optical  thickness  can 
be  varied  from  very  small  to  semi-infinite.  The  single  scattering  albedo 
and  phase  function  give  considerable  insight  into  the  phvsics  governing 
radiative  transfer  within  the  ice.  The  applications  of  the  four  stream 


model  are  discussed  in  Chapter  Five. 


CHAPTER  THREE 


EXPERIMENTAL  APPARATUS  AND  TECHNIQUES 

3.1  APPARATUS  AND  PROCEDURES 

The  experiments  were  performed  in  a  cold  room  whose  air  temperature 
could  be  controlled  over  the  range  -40°C  to  +10°C.  Freezing  occurred  in 
a  cylindrical  plexiglas  tank,  61  cm  high  and  116  cm  in  diameter 
(Figure  3.1).  A  cylindrical  design  was  chosen  to  avoid  stress  concentra¬ 
tions,  enabling  the  use  of  a  very  thin,  .375  cm,  plexiglas  wall  while 
maintaining  structural  integrity.  Mirrored  plexiglas  was  selected  to 
reduce  optical  edge  effects  by  making  the  ice  sheet  appear  larger  in 
areal  extent.  This  was  necessary  to  justify  the  theoretical  assumption 
of  a  plane-parallel  medium.  The  tank's  base  was  1.25  cm  thick  plexiglas 
painted  black  to  minimize  reflection.  To  prevent  leaks  from  the  tank 
silastic  was  used  as  an  inner  sealer.  The  bottom  of  the  tank  was 
insulated  with  7.6  cm  of  "blue  foam".  Insulation  for  the  walls  was 
provided  by  5.1  cm  of  "white  foam"  surrounded  by  5  to  10  cm  of  loose 
zonolite  aggregate  and  topped  by  weatherstripping.  The  tank  and  insula¬ 
tion  were  contained  in  a  wooden  catchbasin  protecting  against  water 
spills  should  a  leak  occur. 

The  source  of  illumination  was  a  "sun"  consisting  of  nine,  500  watt 
incandescent  light  bulbs.  Blue  frosted  superfloods  with  a  color 
temperature  of  4800°K  were  used  since  they  provided  more  short  wavelength 
light  than  conventional  incandescent  bulbs.  To  produce  uniform  illumina¬ 
tion  the  lights  were  mounted  on  a  glossy  white  circular  piece  of  wood, 
one  meter  in  diameter,  in  an  array  consisting  of  a  hexagon  with  an  inner 
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triangle  (Figure  3.2a).  This  configuration  provided  maximum  illumination 
with  minimum  surface  melting.  A  sheet  of  opal  plexiglas  placed  between 
the  "sun"  and  the  ice  surface  acted  as  a  diffuser,  creating  lighting 
conditions  similar  to  cloudy  skies.  An  underwater  lamp  installed  on  the 
bottom  of  the  tank  enabled  the  ice  to  be  viewed  in  transmitted  light. 

The  profiling  spectrophotometer  described  by  Roulet  et  al.  (1974) 
was  used  to  measure  spectral  irradiances  over  the  wavelength  region 
400  to  1000  nanometers.  The  key  feature  of  this  spectrophotometer  was 
the  use  of  fiber  optics  to  transmit  light  from  the  interior  of  the 
medium  to  the  instrument,  allowing  _in  situ  optical  measurements  to  be 
made.  Glass  fiber  optics  probes  were  selected  for  high  transmission 
over  the  entire  wavelength  region  of  interest.  The  fiber  optics  probes 
were  mounted  on  a  central  plexiglas  column  in  a  spiral  array  as  shown  in 
Figure  3.1.  Typically  six  probes  were  used.  An  above  surface  probe  was 
used  for  measurements  of  incident  and  reflected  irradiances,  while  four 
probes  frozen  into  the  ice  at  depths  of  3,  6,  10,  and  19  cm  and  a  sixth 
probe  always  below  the  ice  at  a  depth  of  33  cm,  measured  transmitted 
irradiances.  A  specially  designed  box  (Figure  3.2b)  coupled  the  fiber 
optics  probes  to  the  spectrophotometer.  Tests  were  conducted  to  insure 
that  the  box  allowed  the  probes  to  be  coupled  in  a  reproducible  fashion. 
Optical  data  were  recorded  on  tape  using  an  analog  to  digital  converter. 

To  relate  the  optical  measurements  to  the  physical  state  of  the  ice, 
measurements  of  ice  thickness,  temperature,  and  salinity  were  made. 

Ice  thickness  was  measured  during  the  experiments  using  a  thickness 
gauge  consisting  of  plastic  pipe  (CPVC) ,  .64  cm  in  diameter,  45.7  cm 
long,  and  graduated  in  intervals  of  .25  cm,  with  a  flat  plexiglas  base. 


Frozen  vertically  into  the  ice,  measurements  were  made  by  pulling  the 
gauge  upwards  until  the  base  reached  the  bottom  of  the  ice. 

A  vertical  chain  of  thermistors  was  used  to  obtain  temperature  data. 
Temperatures  were  usually  recorded  at  3  cm  above  the  surface,  at  the 
surface,  and  at  depths  of  2,  4,  6,  8,  10,  12,  14,  16,  18,  20,  22,  and 
34  cm.  The  thermistors  were  potted  in  silastic  on  a  plexlglas  strip  and 
painted  with  marine  varnish  to  provide  electrical  insulation.  Thermal lv 
insulating  effects  of  the  varnish  were  unimportant  since  ice  temperature 
changes  were  slow  compared  to  thermistor  response  time. 

To  determine  salinity  profiles  several  lengths  of  CPVC  pipe,  1.13  cm 
in  diameter,  were  frozen  into  the  ice.  When  a  salinity  profile  was 
needed  a  tube  was  removed  from  the  ice  providing  an  ice  core  which  was 
then  cut  into  2  cm  sections.  After  the  sections  melted  the  salinity  was 
measured  using  an  ENDECO  Type  102  optical  salinometer. 

After  building  and  assembling  the  pieces  of  apparatus  a  series  of 
tests  were  conducted  to  check  the  accuracy  of  the  equipment  and  eliminate 
design  flaws.  The  accuracy  of  the  thickness  gauge  and  the  salinity 
probes  was  checked.  As  a  result  of  these  tests  two  major  additions  to 
the  system  were  made.  A  fan  was  added  to  enhance  growth  by  providing 
air  circulation  and  an  underwater  pump  was  used  to  circulate  the  water, 
thus  retarding  the  growth  of  ice  along  the  wall  of  the  tank.  With  the 
tank  system  it  was  possible  to  grow  ice  25  to  30  cm  thick  similar  to 
that  found  in  the  Arctic,  determine  the  state  and  structure  of  the  ice, 
and  make  in  situ  optical  measurements. 

The  early  testing  of  the  tank  provided  insights  regarding  methods 
of  conducting  successful  experiments  which  were  formalized  into  a 
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standard  set  of  procedures  used  to  perform  the  tank  experiments.  Before 
each  experiment  the  tank  was  thoroughly  cleaned,  the  water  filtered 
through  a  Water  Pic  filter,  and  the  NaCl  and  water  mixed  in  an  inter¬ 
mediate  storage  tank  so  that  any  residue  was  not  added  to  the  tank.  As 
is  frequently  done  in  laboratory  experiments  [Weeks  (1962),  Martin  (1974), 
Eide  and  Martin  (1975),  Davis  and  Munis  (1973),  Lane  (1975),  and  others] 
a  NaCl  solution  was  substituted  for  seawater.  The  cold  temperatures 
during  the  experiments  tended  to  cause  wear  and  breakage  of  the  equip¬ 
ment.  In  particular,  fiber  optics  probes,  thermistors,  the  thickness 
gauge,  and  the  tank's  waterproofing  were  checked  frequently  for  necessary 
repairs  or  replacements.  Due  to  their  short  life,  new  light  bulbs  were 
installed  in  the  "sun"  at  the  beginning  of  every  experiment. 

Each  experiment  consisted  of  two  parts,  a  growth  phase  where  the 
air  temperature  of  the  cold  room  was  held  constant  and  a  warming  phase 
where  the  air  temperature  was  changed  and  the  ice  allowed  to  attain  a 
linear  temperature  profile  at  the  new  air  temperature.  A  group  of 
optical  and  physical  measurements  for  a  given  ice  thickness  and  air 
temperature  is  referred  to  as  a  site.  For  example,  a  typical  experiment 
has  a  growth  phase  consisting  of  eight  sites  for  ice  thicknesses  between 
0  and  26  cm  and  a  warming  phase  with  sites  at  air  temperatures  of  -10°C 
and  0°C.  The  last  site  of  the  growth  phase  is  also  included  in  the 
warming  phase.  Optical  measurements  consist  of  incident,  reflected,  and 
transmitted  irradiances;  while  ice  thickness,  salinity,  temperature,  and 
photographs  constitute  the  physical  measurements.  Table  3.1  outlines 
the  procedure  involved  for  a  particular  site. 


Table  3.1 


Experimental  Procedures 


STEP 

TIME 

(minutes) 

ACTION 

1 

0 

Turn  on  electronics,  warm  up. 

2 

20 

Measure  thermistor  resistances. 

3 

23 

Turn  on  sun. 

4 

30 

Measure  spectrophotometer 
background  signal. 

5 

32 

Make  optical  measurements  of 
incident,  reflected,  and  several 
transmitted  intensities. 

6 

40 

Measure  spectrophotometer 
background  signal. 

7 

41 

Turn  off  sun  and  electronics. 

8 

45 

Measure  thickness. 

9 

50 

Remove  salinity  core,  section, 
and  allow  to  melt  for  later  use. 

10 

60 

Photograph  ice  surface. 

»a 


At  the  end  of  the  growth  phase,  a  core  was  removed  for  photographic 
analysis  of  the  ice  structure.  A  vertical  slab  approximately  two  centi¬ 
meters  thick  was  prepared  and  photographed  in  both  reflected  and  trans¬ 
mitted  light.  A  close-up  profile  was  made  using  transmitted  light.  To 
display  the  crystal  structure  crossed-polaroid  pictures  of  both  hori¬ 
zontal  and  vertical  thin  sections  were  taken  and  the  remainder  of  the 
ice  core  was  stored  for  reference. 

Though  primarily  influenced  by  the  air  temperature  during  freezing, 
the  duration  of  an  experiment  also  depended  on  such  factors  as  the  length 
of  the  warming  phase  and  the  amount  of  breakage  occurring.  Typical  lv, 
experiments  could  be  initiated  at  four  week  intervals. 

3.2  DATA  REDUCTION 
3.2.1  Optical  Measurements 

The  reduction  and  analysis  of  the  optical  information,  as  outlined 
in  Figure  3.3,  was  a  multistep  procedure.  Raw  data,  consisting  of 
spectrophotometer  slit  positions  and  irradiances  were  recorded  on  a 
cassette  tape  recorder.  The  tape  was  read  through  a  software  program  on 
a  Raytheon  704  computer  which  encoded  the  information  generating  a  CDC 
compatible  magnetic  tape.  The  reduction  continued  on  a  CDC  6400 
computer  where  the  data  were  partitioned  into  individual  scans,  consist¬ 
ing  of  about  230  ordered  pairs  of  wavelength  and  irradiance  values  over 
the  wave-length  interval  400  to  1000  nm.  Repeated  scans  for  a  given 
probe  and  site  were  then  averaged  together  and  transposed  onto  a 
standard  wavelength  grid. 
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Since  several  fiber  optics  probes  were  used  concurrently,  it  was 
necessary  to  determine  the  relative  sensitivity  of  each  probe.  Scaling 
factors  Sn  were  selected  for  each  probe  using  data  from  the  no  ice  (open 
water)  measurements  and  the  relationship 


V*) 


F+(n,X)[l-ct(A)] 


(3.1) 


where  F^(n,A)  is  the  transmitted  irradiance  at  probe  n.  Equation  (3.1) 
could  only  be  used  at  a  wavelength  where  no  appreciable  absorption  occurs 
in  the  water.  Smith  (1973),  at  a  wavelength  of  490  nm,  determined  an 
extinction  coefficient  of  .0444  m-1  for  Arctic  Upper  Water.  Substituting 
this  value  into  the  Bouguer-Lambert  model  for  typical  probe  depths 
indicated  that  the  maximum  error  introduced  by  Equation  (3.1)  at  490  nm 
was  less  than  1.5%.  Measured  spectral  irradiances  for  each  probe  were 
then  multiplied  by  the  appropriate  scaling  factor. 

With  scans  on  a  common  grid  and  appropriate  scaling  factors  deter¬ 
mined,  comparisons  were  made  between  different  depths.  Albedos  were 
calculated  from  incident  and  reflected  scans.  Albedos  and  transmitted 
values  were  input  to  the  Dunkle  and  Bevans  photometric  model  for  deter¬ 
mination  of  extinction  coefficients.  Light  transmission  profiles  were 
obtained  from  transmitted  scans  for  various  ice  thicknesses. 

To  aid  in  the  analysis  of  interrelationships  between  optical 
properties  and  other  ice  parameters,  bulk  albedos  defined  by  Equation 
(3.2)  were  calculated  ,-1000 


rx  vj 

J  a  F  (A)dA 

,  ,  4Q0  A  0V _ 

B  .1000 

J  Fn(X)d,\ 

400  0 


(3.2) 
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To  obtain  tlio  necessary  spectral  albedos,  exper  imcnt  a  I  results  ovei  t  ho 
interval  500  to  1000  nm  were  smoothed  using  a  liip.li  order  polynomial  I  it  . 
One  to  the  typically  small  signal  to  noise  ratio  for  wavelengths  loss 
than  500  nm,  experimental  values  were  discarded  and  the  value  at  S0(l  nm 
assumed  to  apply  to  the  entire  400  to  500  nm  region.  This  approximation 
was  based  on  the  results  of  Grenfell  and  Mavkut  (ll>/g)  which  showed  >1.0 
to  be  fairly  constant  below  500  nm.  Values  of  K.,  (V)  for  both  deal  and 
cloudy  skies  were  obtained  from  Sauberer  and  Ournhirn  (ll>58)  and  Cast 
( I960)  and  are  plotted  in  Figure  3.4. 

3.2.2  Physical  Measurements 

Ice  temperatures,  salinities,  and  thicknesses  were  routinely 
measured  during  the  experiments.  Thermistor  readings  ot  resistance  were 
converted  into  temperatures.  Using  temperature  and  salinity  profile 
data  as  input  brine  volumes  were  calculated  from  the  expression 

VH  “  “,055  f  (1.3) 


where  v  is  brine  volume,  S  is  the  salinity  in  parts  per  thousand,  and 
ti 

T  is  the  temperature  in  degrees  Centigrade. 

Ice  thicknesses  were  plotted  as  a  function  of  time  (t  )  and  Ktpi.it  ton 
(3.4)  was  used  to  determine  point  to  point  growth  rates  (f) 
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Theoretical  thicknesses  and  growth  rates  were  calculated  using  a  Stefan 


RELATIVE  IRRADIANCE 


Figure  3.4:  Normalized  Incident  radiation  under:  - clear  skies 

(Sauberer  and  Clrmhlrn,  I9.r)8),  and - cloudy  skies  (Cast,  19t>0). 
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growth  model,  which  had  been  modified  to  use  empirical  formulas  for  the 
heat  of  fusion  and  conductivity  of  sea  ice,  and  compared  to  measured 
values. 


3.3  ERROR  ANALYSIS 
3.3.1  Optical  Measurements 

The  optical  measurements  contained  three  potentially  significant 
sources  of  error:  temporal  and  spatial  variation  in  illumination, 
instrumental  noise,  and  the  bandpass  of  the  instrument.  In  designing 
and  testing  the  experimental  conf iguration  several  steps  were  taken  to 
minimize  errors  caused  by  the  system  of  illumination.  The  sun  was 
enlarged  from  six  bulbs  to  nine,  thus  providing  an  incident  radiation 
field  of  greater  areal  uniformity.  In  early  tank  testing  the  incident 
radiation  was  measured  at  the  beginning  and  end  of  each  set  of  readings 
to  examine  the  variability  of  the  lighting.  It  was  determined  that 
there  was  no  appreciable  change  in  incident  irradiance  during  a  set  of 
measurements,  but  that  the  variation  between  sites  was  large  enough  to 
necessitate  an  incident  measurement  for  each  site. 

The  electronic  noise  level  was  independent  of  wavelength  and  equal 
to  1  2.5  mV  for  the  more  sensitive  instrument  setting  and  1  1  mV  for  the 
less  sensitive.  However,  since  this  noise  level  affected  a  signal  which 
varied  in  magnitude  with  wavelength,  the  percent  error  had  a  decided 
wavelength  dependence.  Figure  3,5  illustrates  the  resultant  ranges  in 
calculated  albedos  for  representative  cases.  Due  to  low  signal  levels 
the  uncertainty  between  400  and  450  nm  was  so  extreme  as  to  preclude  any 
meaningful  interpretation  of  the  data. 


ALBEDO 
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Figure  3.5:  Range  of  uncertainty  in  albedo  due  to  instrumental  noise: 
The  lower  curve  is  ice  11.7  cm  thick  at  an  air  temperature  of  -10°C 
and  the  upper  curve  is  ice  14.6  cm  thick  at  an  air  temperature  of 
-30°C. 
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Spectral  errors  were  Introduced  at  long  wavelengths  by  the  large 
bandpass  of  the  instrument.  This  bandpass  error  was  manifested  as  an 
apparent  decrease  in  the  wavelength  dependence  of  a  property.  The 
spectrophotometer,  using  a  glass  prism,  had  less  resolving  power  in  the 
red  portion  of  the  spectrum  than  in  the  blue.  Also,  as  was  shown  in 
Figure  1.1,  previous  research  indicated  an  accelerating  increase  in 
extinction  coefficient  with  wavelength  for  sea  ice.  Thus,  the  loss  of 
spectral  resolution,  while  not  important  at  shorter  wavelengths,  became 
critical  at  longer.  Deconvolution  techniques  were  used  in  an  attempt  to 
increase  spectral  resolution.  However,  deconvolution  magnified  small 
signal  variations  into  large  oscillations  (Figure  3.6)  yielding  unrealis¬ 
tic  values  and  making  it  inapplicable  in  this  experiment.  Since  its 
spectral  dependence  was  weaker,  albedo  was  not  as  strongly  influenced 
by  this  error  as  extinction  coefficient.  Because  of  the  bandpass  broad¬ 
ening,  extinction  coefficients  calculated  for  wavelengths  greater  than 
800  nm  were  not  reported. 

3.3.2  Physical  Measurements 

Thicknesses  and  the  depths  of  the  fiber  optics  probes  were  measured 
with  an  accuracy  of  *  1.5  mm.  There  was  no  preferential  growth  of  ice 
along  the  thickness  gauge.  Determining  salinity  profiles  the  error  in 
reading  the  optical  sallnometer  was  1  .1  °/00  and  the  sectioning  was 
accurate  to  *  1  mm  in  position.  Additional  uncertainty  in  salinity  was 
caused  by  the  small  sample  volume.  Because  of  brine  drainage,  salinity 
measurements  from  the  bottom  section  of  the  core  were  erroneously  small. 
Thermistor  positions  were  known  to  1  1  mm  while  the  error  in  temperature 
varied  from  .2  to  .  4°C  as  illustrated  in  Figure  3.7. 
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Figure  3.7:  Thermistor  uncertainty. 


CHAPTER  FOUR 


EXPERIMENTAL  RESULTS 


4.1  INTRODUCTION 

A  series  of  18  experiments  were  performed  using  the  apparatus  and 
methods  discussed  in  Chapter  Three.  Six  of  these  occurred  during  the 
testing  phase,  yielding  primarily  information  regarding  improvements  to 
the  experimental  system.  After  the  testing  was  completed,  it  was 
determined  that  eight  cases  would  best  satisfy  the  experimental 
objectives.  Due  to  assorted  electronic  failures,  it  was  necessary  to 
perform  twelve  experiments  to  completely  fulfill  these  data  requirements. 

Experiments  were  conducted  for  different  water  salinities  and  over 
a  range  of  air  temperatures.  On  the  basis  of  water  salinity,  experiments 
were  grouped  into  three  categories;  fresh  water  (Experiments  1  and  2), 
salinity  (31  °/00)  comparable  to  Arctic  Ocean  water  (3-6),  and  inter¬ 
mediate  (16  °/00)  salinity  water  (7,  8).  Salt  ice  was  grown  at  air 
temperatures  of  -10,  -20,  -30,  and  -37°C.  The  environmental  conditions 
for  the  experiments  are  summarized  in  Table  4.1  with  more  detailed 
information  regarding  individual  sites  displayed  in  Table  4.2. 


4.2  SELECTED  RESULTS 
4.2.1  Fresh  Ice  Results 

During  the  first  two  experiments  fresh  ice,  grown  at  an  air  tempera¬ 
ture  of  -15°C,  was  investigated  to  compare  to  the  salt  ice  cases  and  to 
examine  the  influence  of  air  bubbles.  Experiments  1  and  2  exhibited  a 
range  of  bubble  conditions  including  bubble  free  ice,  spherical  bubbles, 
and  long  cylindrical  vertically  oriented  bubbles. 
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Table  4.1 

Summary  of  Tank  Experiments 


Experiment 

Number 

Water 

Salinity 

(O/  ) 

i  '  oo> 

Air  Temperature 
During  Freezing 
(°C) 

Probe  Depths 
(cm) 

Number  of  Sites 
Growth:  Warming 

1 

0 

-15 

9.7,  17.7,  18.2 
18.2,  33.0 

!  8:0 

2 

0 

-15 

3.0,  10.7,  19. C 
34.0 

)  5:0 

3 

30.6 

-10 

3.1,  6.9,  11.0 
19.3,  34.5 

7:2 

4 

30.2 

-20 

2.5,  6.3,  10.5 

18.5,  33.6 

7:3 

5 

30.8 

-30 

3.2,  6.9,  11.0 
19.4,  34.4 

6:2 

6 

31.1 

-37 

3.0,  6.4,  10.9 
19.0,  34.1 

7:2 

i 

7 

16.3 

-10 

3.2,  6.7,  11.0 
19.9,  34.5 

7:3 

8 

16.5 

-20 

3.0,  6.4,  10.9, 
20.3,  34.3 

,  7:2 
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Table  4.2 

Ice  Thickness  (cm)  and  Air  Temperature  (°C) 
for  Experimental  Sites 


Experiment 


1 

2 

3 

4 

5 

6 

7 

8 

Open 

Open 

Open 

Open 

Open 

Open 

Open 

Open 

Water 

Water 

Water 

Water 

Water 

Water 

Water 

Water 

2.54 

9.2 

4.3 

2.0 

3.2 

3.8 

4.5 

3.6 

-15 

-15 

-10 

-20 

-30 

-32 

-20 

-10 

4.4 

12.5 

7.6 

4.6 

11.2 

7.6 

9.3 

7.9 

-15 

-15 

-10 

-20 

-30 

-32 

-20 

-10 

6.4 

18.2 

11.7 

7.6 

14.6 

14.5 

12.1 

11.6 

-15 

-15 

-10 

-20 

-30 

-33 

-20 

-10 

8.9 

22.4 

18.5 

11.7 

19.6 

17.5 

15.4 

14.2 

-15 

-15 

-10 

-20 

-30 

-36 

-20 

-10 

10.4 

— 

21.1 

17.8 

23.9 

23.5 

19.3 

19.6 

-15 

-10 

-20 

-30 

-37 

-20 

-10 

12.2 

— 

27.2 

21.6 

29.0 

26.0 

22.6 

23.0 

-15 

-10 

-20 

0 

-37 

-20 

-10 

16.8 

— 

28.4 

25.7 

29.2 

28.0 

23.8 

23.4 

-15 

0 

-16 

-20 

-22 

0 

0 

- 

- 

25.9 

-5 

26.2 

-10 

- 

29.0 

_2 

23.8 

-20 

23.4 

-10 

- 

— 

- 

27.9 

0 

- 

- 

- 

23.4 

-20 
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The  top  seven  centimeters  of  ice  grown  in  Experiment  1  were 
virtually  bubble  free.  Bubbles  then  appeared  and  gradually  increased 
in  number  between  7  and  12  cm,  with  a  high  density  of  bubbles  being 
present  from  12  cm  to  the  bottom  of  the  ice.  The  bubbles  were  predomi¬ 
nately  spherical  with  a  few  being  cylindrical. 

Measurements  across  the  bubble  free  region  indicated  that  within 
the  limits  of  error  bubble  free  ice  and  water  are  optically  indistin¬ 
guishable.  This  is  in  agreement  with  other  reported  measurements  of 
extinction  coefficient  for  pure  ice  and  water  (Goodrich,  1970).  Based 
on  the  results  of  previous  researchers,  transmission  differences,  for  7 
cm  thick  ice,  between  the  two  cases  are  calculated  to  be  no  greater 
than  1.5*  for  wavelengths  between  400  and  800  nm;  a  difference  less  than 
the  experimental  error  levels. 

In  the  bubbly  region  of  Experiment  1  albedo  evidenced  a  slight 
increase  with  thickness.  Figure  4.1  illustrates  spectral  albedos 
measured  for  fresh  ice  of  thickness  6.4,  12.2,  and  16.8  cm.  While  the 
albedo  of  bubble  free  ice  (Curve  1)  is  the  same  as  that  for  water,  there 
is  an  increase  as  the  ice  thickens  (Curves  2  and  3).  The  increasing 
density  of  air  bubbles  in  the  ice  results  in  greater  scattering  and  a 
larger  albedo.  There  was  no  noticeable  difference  in  extinction 
coefficient  between  these  cases.  However,  the  ability  to  measure 
variations  in  extinction  coefficient  is  less  precise  due  to  additional 
errors  introduced  by  the  inter-probe  scaling  factors. 

Experiment  2  exhibited  a  very  different  type  of  bubble  structure 
with  vertical  strings  of  small  bubbles  and  long,  thin,  vertically 
oriented  cylindrical  bubbles  prevalent.  The  diameter  of  the  cylinders 
was  roughly  .2  mm  with  lengths  ranging  from  2  to  40  mm;  7  mm  being 
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Figure  4.1:  Spectral  albedos  for  fresh  ice  of  thickness: 
(1)  6.4  cm;  (2)  12.2  cm;  and  (3)  16.8  cm. 


typical.  In  many  instances  closer  inspection  revealed  that  what 
appeared  to  be  one  cylindrical  bubble  was  a  linear  array  of  many  small 
bubbles.  However,  in  a  general  sense,  tt.«  hubbies  can  be  described  as 
long,  narrow,  vertically  orient  d  cvi-ttder*.  Spectrophotometer  results, 
even  for  ice  22.4  cm  thick,  indicated  n.  dtfteteu. <  in  either  albedo 
or  extinction  coefficient  from  t  of  ,t 

Qualitatively,  however,  a  vers  Inti  •  t  n  •  •  .»»  ure  was  observed. 

The  photographs  in  Figure  4.2,  .  •:  ‘  ■  .  .  t  .  of  illumination, 

were  taken  under  identical  cond.ti.  s.  ; .  .to  .. .  ;  t.\ )  the  ice  is 
illuminated  from  above,  parallel  \  'd«-r  .  .here  appears  to 

be  little  scattering.  When  illuminated  t t or  ,d. ,  perpendicular  to 

the  cylinders,  as  in  Figure  a. 2  lb  1 ,  the  increase  ta  scattering  is 
readily  apparent.  The  amount  of  scattering  due  to  an  air  bubble  can  be 
expected  to  be  roughly  proportional  to  the  projected  area  of  the  bubble 
onto  a  plane  perpendicular  to  the  direction  of  illumination.  A  cylindri¬ 
cal  air  bubble  has  a  different  projected  area  perpendicular  to  and 
parallel  to  its  axis  of  orientation.  For  incident  light  parallel  to  the 
cylinders'  orientation  the  projected  area  of  the  bubble  is  nr  and  is 
2fr  for  light  perpendicular  where  t  is  the  cylinder  length  and  r  is  the 
cylinder  radius.  The  ratio  of  the  projected  areas  is 


For  r  equal  to  . 1  mm  and  i  equal  to  7  mm,  typical  values  for  this 
experiment,  is  .0225,  implying  that  the  scattering  should  be  on 

the  order  of  45  times  greater  for  the  perpendicular  case.  Further  quan¬ 
titative  Investigations  of  this  phenomena  are  suggested  in  Chapter  Six. 
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g  by  cylindrical  bubbles  in  fresh  ice: 
parallel  to  cylinders,  and  (B)  illumination 


A. 2. 2  Salt  Ice  Results:  Experiment  3 

Experiments  3  through  6  are  designed  to  determine  the  optical 
properties  of  young  thin  salt  ice  under  a  range  of  conditions.  To  avoid 
a  tedious  itemizing  of  data,  only  the  results  from  one  representative 
experiment.  Experiment  5,  are  discussed  in  detail.  Results  from  the 
remaining  experiments  are  summarized  in  Appendix  A. 

The  initial  water  salinity  in  Experiment  3  was  30.8  During 

the  growth  phase,  the  air  temperature  was  held  at  -30°C.  Measurements 
were  made  at  ice  thicknesses  of  0.0,  3.2,  7.6,  14.6,  19.6,  and  23.‘>  cm. 
After  completing  the  growth  phase  the  air  temperature  was  raised  to  -2°C 
and  measurements  were  taken  after  allowing  sufficient  time  for  the  Ice  to 
achieve  a  linear  temperature  gradient.  For  the  final  site  the  air 
temperature  was  lowered  to  -20°C  and  the  process  repeated.  Figures  4.3 
and  4.4  display  temperature  and  salinity  profiles  recorded  in  conjunction 
with  the  optical  measurements.  Temperature  profiles,  within  thermistor 
error  limits,  were  linear.  The  ice  surface  temperature  was  always 
warmer  than  the  air  temperature  of  -30°C,  though  approached  it  as  the 
ice  grew  thicker. 

The  ice  salinities  were  quite  high  with  the  profiles  exhibiting 
considerable  scatter.  The  profiles,  for  the  most  part,  did  not  exhibit 
the  strongly  c-shaped  curve  associated  with  field  measurements  of  young 
sea  ice.  Since  the  amount  of  brine  entrapped  increases  witli  growth 
rate,  the  salinity  of  young  ice  decreases  from  the  surface  downward 
until  reaching  the  skeletal  growth  layer  at  the  bottom  where  due  to  the 
porous  structure,  the  salinity  again  increases.  However,  in  the  experi¬ 
ments  the  salinity  decrease  over  the  middle  region  was  not  as  pronounced 
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4.4:  Salinity  profiles  for  Experiment  5.  The  numbers  correspond  t 
same  cases  discussed  in  Figure  4.3.  The  scale  Is  the  same  for  each 
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since  the  decrease  in  growth  rate  was  counterbalanced  by  an  increase  in 
water  salinity. 

Ice  thickness  as  a  function  of  time  and  the  accompanying  growth 
rates  are  displayed  in  Figure  4.5.  The  growth  r <te  smoothly  decreased 
as  the  ice  thickened.  The  decrease  in  growth  rate  caused  by  the 
increase  in  water  salinity  was  less  than  three  percent. 

At  the  conclusion  of  ttie  growth  phase  a  vertical  cross  section  was 
prepared  and  photographed  in  transmitted  natural  light  (Figure  4.6)  and 
polarized  light  (Figure  4.7).  Viewed  in  natural  light  the  ice  appears 
fairly  homogeneous  with  very  faint  bands  at  depths  of  9,  11,  and  18  cm. 
However,  crossed  polaroids  reveal  a  more  complex  crystal  structure.  At 
the  surface  is  a  very  thin  (.3  cm)  layer  containing  ice  crystals  with 
horizontal  c-axes.  Beneath  that  is  a  zone,  roughly  4  cm  thick,  consist¬ 
ing  of  a  jumble  of  crystals,  where  the  c-axis  gradually  changes  from  a 
vertical  to  a  horizontal  orientation.  From  4  to  8  cm  the  crystals  are 
long  and  ordered  with  horizontal  c-axes.  At  a  depth  of  8  cm  these 
crystals  are  terminated  by  another  zone  of  small  disordered  crystals. 

At  11  cm  the  c-axis  orientation  again  becomes  horizontal  and  remains  so 
throughout  the  rest  of  the  sample. 

Spectral  albedos  for  various  ice  thicknesses  are  plotted  in  Figure 
4.8.  Albedos  from  400  to  450  run  (dotted  lines)  are  assumed  constant. 
There  is  a  pronounced  increase  in  albedo  as  the  ice  becomes  thicker. 
Plotted  for  comparison  (dashed  curve)  are  the  albedos  of  thick,  melting, 
multiyear  white  ice  and  melting,  first-year  blue  ice  (dots  and 
dashes)  taken  from  Grenfell  and  Maykut  (1978).  The  albedo  of  the  multi¬ 
year  ice  is  similar  to  that  for  cold,  growing  ice  only  23.9  cm  thick. 

The  increase  in  albedo  with  ice  thickness  is  wavelength  dependent,  with 
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Ice  growth  data  for  Experiment  5. 


Kiguro  A.  7:  let*  slab  from  Kxpor  imont  ri  phot  ograptu'il  through  urossuil 
po 1 aro ids . 
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Figure  4.8:  Spectral  albedos  for  Experiment  5:  (1)  open  water,  cloudy 
sky  bulk  albedo  of  .07;  (2)  Ice  thickness  of  3.2  cm,  .09;  (3)  7.6  cm 

.31;  (4)  14.6  cm,  .44;  (5)  19.6  cm,  .53;  (6)  23.9  cm,  .63;  - 

melting  first-year  blue  ice  (Grenfell  and  Maykut,  1978);  and  - 

melting  multiyear  white  ice  (Grenfell  and  Maykut,  1978). 
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a  greater  increase  occurring  at  shorter  wavelengths.  Table  4.3  corrob¬ 
orates  this  result,  showing  an  increase  in  the  ratio  O5oo/a800  from 
1.07  at  3.2  cm  to  1.55  at  23.9  cm. 

Table  4.3 

Ratio  of  «5oo/a800  as  a  Function  of  Ice  Thickness 


H 

“500 

a800 

a50C)/a800 

- 

3.2 

.094 

.088 

1.07 

7.6 

.35 

.32 

1.10 

14.6 

.51 

.36 

1.42 

19.6 

.62 

.41 

1.51 

23.9 

.76 

.49 

1.55 

Light  attenuation  within  the  ice  is  both  depth  and  wavelength 
dependent.  This  is  illustrated  by  profiles  of  light  transmission 
(Figure  4.9),  measured  at  the  conclusion  of  the  growth  phase.  The 
attenuation  is  decidedly  nonexponential  with  depth.  Attenuation  is 
greatest  near  the  surface  and  decreases  with  depth.  This  is  probably 
due  to  depth  dependent  structural  variations  in  the  ice  and  is  examined 
in  greater  detail  in  Chapter  5.  Attenuation  becomes  greater  as  wave¬ 
length  increases.  Through  20  cm  ice  20%  of  500  nm  light  is  transmitted, 
while  at  800  nm  only  12%  penetrates. 

Incident,  reflected,  and  transmitted  spectral  intensities  were  used 
in  conjunction  with  the  Dunkle  and  Bevans  model  to  calculate  spectral 
extinction  coefficients.  Values  plotted  in  Figure  4.10  are  for  ice  23.9 
cm  thick,  using  measurements  from  a  depth  of  19.4  cm.  Extinction 


NORMALIZED  TRANSMITTED  RADIANCE 


Figure  4.9:  Spectral  transmission  profiles  from  Experiment  5  for 
23.9  cm  thick  ice  with  an  air  temperature  of  -30°C:  (1)  500  nm; 

(2)  600  nm;  (3)  70 0  nm;  and  (4)  800  nm. 
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exhibits  a  marked  wavelength  dependence,  increasing  slowly  between  450 
and  600  nm,  then  sharply  from  600  to  800  nm.  At  600  nra  the  value  is 
only  about  20%  larger  than  at  500  nm,  but  by  800  nm  it  is  90%  larger 
than  the  value  at  500  nm. 

Extinction  coefficients  at  all  wavelengths  are  considerably  larger 
than  previously  reported  results  for  thick  ice.  At  500  nm  the  extinction 
coefficient  for  this  case  is  twice  that  measured  in  the  surface  granular 
layer  of  multiyear  white  ice  and  ten  times  greater  than  values  for  the 
interior  of  multiyear  white  ice  (Grenfell  and  Maykut,  1978).  The  bulk 
extinction  coefficient,  calculated  over  the  wavelength  region  400  to  800 
nm,  is  .084  cm-1.  This  optical  bulk  extinction  coefficient  is  approxi¬ 
mately  five  times  larger  than  the  values  of  .017  and  .012  cm-1  measured 
by  Untersteiner  (1961)  and  Weller  and  Schwertfeger  (1967)  for  multiyear 
ice  and  first-year  ice,  respectively. 

The  reason  for  the  above  differences  must  be  related  to  other 
physical  properties  of  the  ice.  In  these  experiments  the  ice  is  thinner, 
colder,  and  saltier  than  the  sea  ice  observed  in  previous  work.  The 
causes  of  larger  extinction  coefficients  in  the  thin  cold  ice  case  are 
examined  in  Section  4.3. 

4.2.3  Thermal  Determination  of  Extinction  Coefficients 

In  three  of  the  experiments  an  independent  analysis  of  attenuation 
in  the  ice  is  undertaken  to  provide  supplementary  confirmation  of  the 
extinction  coefficients.  A  method  relating  the  heat  absorbed  in  the  ice 
to  the  extinction  coefficient,  comparable  to  that  of  Untersteiner  (1961), 


is  used. 
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The  ice  is  illuminated  for  200  minutes,  with  temperature  profiles 
taken  at  the  start  and  finish  of  that  period.  Knowing  these  temperature 
profiles  and  the  salinity  profile,  the  absorbed  heat  as  a  function  of 
depth  is  calculated.  From  basic  thermodynamics 

dQ  =  cgpdT  (4.2) 

where  dQ  is  the  heat  absorbed  per  unit  volume,  c  =  0.5  +  4.1(S/T‘>) 

s 

(Untersteiner ,  1961),  and  p  is  the  density  =  p.  (l-vD)  +  p  v  *  p 

ICG  D  W  D  ICG 

Integrating  Equation  (4.2)  over  temperature, 

Q2-Qj  -  p(T2-T1)(0.5  +  4.1  -)  .  (4.3) 

Heat  is  transmitted  within  the  ice  by  both  conduction  and  absorption 
of  radiation.  A  computer  program,  using  a  forward  differences  treatment 
of  the  heat  conduction  equation  (Richtmyer,  1957),  was  written  to 
evaluate  the  contribution  of  conduction  to  temperature  changes  within 
the  ice.  For  the  experimental  conditions  conduction  was  found  to  be 
insignificant  below  a  depth  of  four  centimeters,  thus  heating  in  this 
region  is  due  only  to  absorbed  radiation.  Consequently, 

3F 

AQi  =  a^;At  »  (A-4) 

for  time  interval  At  at  some  depth  z^.  Since  attenuation  over  a  small 
z  interval  can  be  approximated  by  a  Bouguer-Lambert  model,  Equation  (4.4), 
at  depths  z^  and  is  expressed  as 


AQ4  =  -  kFqc"*2*  At 


(4.5) 
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AQi+i  =  -  KF0e'K2i+l  At  . 
Dividing  Equation  (4.5)  by  (4.6)  and  solving  for  k. 


K 


AQj+i 

AQi 


/ 


(zi+r2i) 


(A. 6) 


(A. 7) 


The  k  calculated  by  Equation  (A. 7)  is  a  bulk  extinction  coefficient, 
integrated  over  the  spectrum  available  at  depth  z .  Figure  A, 11  illu¬ 
strates  salinity  and  temperature  profiles,  before  and  after  radiative 
heating  for  ice  26  cm  thick  and  grown  at  an  air  temperature  of  -37°C. 

The  heat  absorbed  in  two  centimeter  intervals  for  depths  from  0  to  20  cm 
was  calculated  and  displayed  in  Figure  A. 12.  For  greater  accuracy,  a 
best  fit  of  the  curve  was  determined  for  the  region  from  A  to  12  cm. 
Shallower  depths  were  influenced  by  conduction,  while  small  temperature 
changes  at  larger  depths  increased  errors.  The  bulk  extinction  coeffi¬ 
cient  determined  in  this  fashion  was  0.19  cm-1,  in  general  agreement 
with  the  optical  results. 


A. 3  CORRELATION  OF  OPTICAL  PROPERTIES  AND  PHYSICAL  PARAMETERS 
A . 3 . 1  Overview 

Sea  ice  is  a  complex  medium  with  a  state  and  structure  dependent 
on  such  parameters  as  thickness,  growth  history,  temperature,  salinity, 
and  brine  volume.  One  of  the  goals  of  this  research  is  to  investigate 
correlations  between  the  optical  properties  and  the  physical  state  and 
structure  of  the  ice.  Complicating  this  task  are  the  interrelationships 
among  the  physical  parameters  themselves;  for  example,  salinity  and 
brine  volume  are  explicitly  related  by  Equation  (3.2).  More  difficult 


Klfiurc  4.11:  IVmporai  ut  o  ami  salinity  pro!  lies  useil  In  thermal 
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Figure  4.12:  Profile  of  heat  absorbed  during  a  200  minute 
period  of  illumination. 
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to  define  is  the  dependence  of  salinity  on  growth  rate  and  changes  in 
temperature  profile  as  the  ice  grows.  The  analysis  is  further  compli¬ 
cated  by  the  variability  of  ice  surface  conditions. 

To  quantitatively  determine  the  influence  of  an  individual  para¬ 
meter  it  is  necessary  either  to  perform  experiments  where  effects  of  the 
parameter  are  isolated  or  to  apply  multivariate  analysis  techniques  to 
the  data.  Unfortunately,  due  to  the  nature  of  the  experiments,  it  is 
impossible  to  strictly  isolate  a  single  variable.  Insufficient  data 
prohibit  the  meaningful  application  of  multivariate  analysis.  However, 
by  comparing  a  number  of  cases  under  varying  conditions  a  qualitative 
understanding  of  the  correlation  between  the  optical  properties  and  the 
physical  parameters  can  be  achieved. 

4.3.2  The  Effects  of  Ice  Thickness 

For  increasing  ice  thickness  the  albedo  increases  and  the  trans¬ 
mission  decreases.  As  was  indicated  earlier,  the  influence  of  thickness 
on  the  optical  properties  may  be  in  part  due  to  changes  in  the  ice 
temperature  profile.  Figure  4.13  displays  isopleths  of  albedo  as  a 
function  of  wavelength  and  thickness  from  the  growth  phase  of  Experiment 
4  =  -20°C).  As  in  the  other  experiments,  the  increase  of  albedo 

with  thickness  was  observed  to  be  greatest  at  short  wavelengths. 

The  thickness  dependence  of  bulk  albedo  for  the  growth  phases  of 
Experiments  3-6  is  illustrated  in  Figure  4.14;  for  comparison  bulk 
albedos  measured  by  Weller  (1972)  using  a  Kipp  radiometer  (400  to  3000  nm) 
are  also  Included.  The  bulk  albedo  curves  exhibit  an  initial  rapid 
increase  followed  by  a  gradual  tapering  as  the  thickness  increases .  The 
large  increase  in  Curve  4,  as  the  thickness  increases  from  17.4  to  23.7 


BULK  ALBEDO 


ICE  THICKNESS  (cm) 


Figure  4. 14:  Bulk  albedo  as  a  function  of  ice  thickness  for: 

(1)  Experiment  3  (air  temperature  of  -10°C) ;  (2)  Experiment  A  (-20°C) ; 

(3)  Experiment  5  (-30°C);  (A)  Experiment  6  (-37°C);  - - - - - 

Weller,  1972;  and  - theoretical  curve  using  Equation  5.1. 
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cm,  is  due  Co  the  temperature  in  the  upper  few  centimeters  of  the  ice 
dropping  below  the  eutectic  point,  resulting  in  the  precipitation  of 
solid  salts. 

Using  a  simple  approach  a  relationship  between  bulk  albedo  and 
thickness  is  derived.  Assuming  single  scattering,  vertically  homo¬ 
geneous  reflectivity,  exponential  attenuation  within  the  ice,  and 
specular  reflection  of  magnitude  RQ  ,  the  contribution  to  the  albedo  from 
a  layer  a  deep  and  dz  thick  is 

a  (dz)  =  be“2az  .  (4.8) 

D 

Integrating  over  z  from  0  to  H  gives 

- 

“b(H)  -  h  (1_e_2aH)  +  R0  (4  .°) 

where  the  specular  reflection  Rp  is  seven  percent,  b  and  a  are  the 
amounts  backscattered  and  absorbed  in  a  layer  dz.  Thickness  dependent 
bulk  albedos  from  Experiment  3  conducted  at  an  air  temperature  of  -10°C 
are  least  influenced  by  changes  in  temperature  profiles  and  are  selected 

| 

to  fit  Equation  (4.9).  A  nonlinear  regression  program,  using  Equation 
(4.9)  and  ordered  pairs  of  bulk  albedo  and  thickness  as  input,  calculated 
b  and  a  to  be  .0184  and  .0149,  respectively,  with  a  correlation  coeffi¬ 
cient  of  .98.  Not  only  does  Equation  (4.9)  fit  the  data  reasonably  well 
(dashed  curve,  Figure  4.14),  but  in  the  limiting  case  of  very  large  ice 
thickness  (H**)  it  gives  a  physically  realistic  value  of  .69. 

A  relationship  such  as  Equation  (4.9),  while  offering  little 
physical  insight,  is  of  practical  value.  Maykut  (1977)  discusses  the 
importance  of  thickness  variation  in  the  large-scale  (regional)  heat  and 
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mass  balance  of  the  Lee  pack.  If  the  ice  thickness  distribution  g(H)  is 
known,  then  the  regional  average  of  some  thickness  dependent  property 
4>(H)  is  given  by 

/  max 

+  =  I  <KH)g(H)dH  .  (A. 10) 

H  , 
min 

Substituting  Equation  (4.9)  into  Equation  (4.10)  for  a  known  thickness 
distribution  of  bare  ice,  a  regional  albedo  average  can  be  calculated. 


4.3.3  The  Effects  of  Air  Temperature  During  Freezing  and  Growth  Rate 
For  Experiments  3-6  the  air  temperatures  during  the  growth  phase 
were  -10,  -20,  -30,  and  -37°C,  enabling  an  investigation  of  how  differ¬ 
ences  in  the  air  temperature  during  freezing  and  thus  the  growth  rate 
affect  the  optical  properties.  According  to  the  Stefan  equation  (Carslaw 
and  Jaeger,  1959)  the  air  temperature  and  growth  rate,  assuming  a  linear 
temperature  gradient  in  the  ice,  are  related  by 


f 


dH 

dt 


Vice11 


(4.11) 


where  is  the  latent  heat  of  fusion,  Tfl  is  the  air  temperature  (°K) , 

Tj  is  the  freezing  temperature,  ct  is  the  coefficient  of  thermal 

conductivity,  p.  is  the  density  of  ice,  and  H  is  the  ice  thickness, 
ice 

Because  of  the  dependence  of  growth  rate  on  ice  thickness,  to  facilitate 
comparisons,  each  experiment  was  represented  by  a  growth  rate  determined 
when  the  ice  was  10  cm  thick. 

Data  from  the  growth  phase  of  Experiments  3-6  were  used  to  construct 
contours  of  bulk  albedo  as  a  function  of  ice  thickness  and  air  tempera¬ 
ture  (Figure  4.15)  and  as  a  function  of  ice  thickness  and  growth  rate 
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Figure  4.15:  Isopleths  of  bulk  albedo  as  a  function  of 
air  temperature  during  freezing  and  ice  thickness. 


THICKNESS  (cm) 


Figure  4.16:  Isopleths  of  bulk  albedo  as  a  function  of 
growth  rate  and  ice  thickness 


(Figure  A. 16).  Lack  of  data  makes  the  portion  of  the  0.5  contour  above 
17  cm  uncertain.  Both  figures  display  the  same  general  characteristic: 
a  weak  dependence  of  albedo  on  air  temperature  and  growth  rate  when  the 
ice  is  thin  which  becomes  increasingly  strong  as  the  ice  thickens.  The 
very  large  albedos  (upper  right  hand  corner  of  plots)  are  a  result  of 
the  formation  of  highly  reflective  solid  salts  when  the  ice  temperature 
dropped  below  the  eutectic  point. 

Spectral  albedos,  selected  from  Experiments  3-6,  for  ice  of 
comparable  thickness  (about  20  cm)  are  plotted  in  Figure  4.17.  In  all 
cases,  the  albedo  increases  as  the  growth  rate  increases  (air  temperature 
decreases),  with  the  albedo  more  than  doubling  between  Curves  1  and  4. 
Displayed  in  Figure  4.18  are  corresponding  spectral  extinction  coeffi¬ 
cients  calculated  for  the  same  cases  shown  in  Figure  4.17.  Ice  grown 
under  colder  conditions  (faster  growth)  exhibits  larger  extinction 
coefficients,  increasing  at  600  nm  by  a  factor  of  five  between  Curves  1 
and  4. 

The  experiments  discussed  so  far,  in  addition  to  different  growth 
conditions,  also  have  dissimilar  ice  temperature  profiles.  Therefore, 
observed  differences  may  be  due  to  variations  both  in  growth  rate  and  in 
temperature  profile.  To  more  fully  isolate  the  effect  of  the  growth 
rates  on  the  optical  properties,  sites  are  compared  where  the  tempera¬ 
ture  profiles  are  nearly  the  same.  In  Experiments  3,  5,  and  6,  as  part 
of  the  warming  phase,  the  ice  was  allowed  to  come  to  thermal  equilibrium 
at  a  temperature  near  the  melting  point.  Though  less  pronounced  than 
results  from  growth  phase,  warming  phase  results  show  a  distinct  correla¬ 
tion  between  growth  rate  and  optical  properties  of  the  ice  (Figures  4.19 
and  4.20).  Table  4.4  compares  relative  and  absolute  values  of  albedo 
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Figure  4.17:  Spectral  albedo  comparisons  for  sea  ice  grown  at  different 
air  temperatures.  The  growth  rates  were  determined  when  the  ice  was 
10  cm  thick. 


Experiment 

Thickness 

(cm) 

Air  Temp 
<°C) 

Growth  Rate 
(x  10-r'cm/sec) 

Bulk  Albedo 
(cloudy) 

(1)  3 

21.1 

-10 

4.5 

.37 

(2)  4 

21.6 

-20 

8.0 

.42 

(3)  5 

19.6 

-30 

11.4 

.53 

(4)  6 

23.5 

-37 

17.0 

.86 

EXTINCTION  COEFFICIENT  (cm”' 


WAVELENGTH  (nm) 


Figure  4.19:  Comparison  of  spectral  albedos  for  ice  of  similar 
temperature  profiles,  but  grown  at  different  air  temperatures. 

In  each  case  the  ice  is  roughly  28  cm  thick  and  the  air 
temperature  is  approximately  0°C:  (1)  Experiment  6,  grown  at  an 

air  temperature  of  -37°C,  and  cloudy  sky  bulk  albedo  of  .51; 

(2)  Experiment  5,  -30°C,  .42;  and  (3)  Experiment  3,  -10°C,  .33. 
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Figure  4.20:  Comparison  of  spectral  extinction  coefficients 
for  ice  of  similar  temperature  profiles,  but  grown  at 
different  air  temperatures.  The  conditions  and  numbers 
correspond  to  the  same  cases  described  in  Figure  4.19. 
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and  extinction  coefficient  obtained  from  Figures  A. 17  through  4.20. 

Even  when  brought  to  the  same  temperature,  faster  grown  sea  ice  exhibits 
larger  albedos  and  extinction  coefficients  than  slower  grown  ice. 


Table  4.4 

Influence  of  Growth  Rate  on  otgoo  and  <*600 


Growth  Phase 
Measurements 

Ice  of  Similar 
Temperature  Profile 

Experiment  Number 

3 

4 

5 

6 

3 

5 

6 

Growth  Rate  at  10  cm 
(x  10" 5  cm/sec) 

4.5 

8.0 

11.4 

17.0 

4.5 

11.4 

17.0 

<*600 

.40 

.49 

.55 

.92 

.36 

.45 

.53 

Normalized  to  Experiment 

3 

1.00 

1.23 

1.38 

2.30 

1.00 

1.25 

1.47 

i<600(cm-1) 

.025 

.031 

.080 

.126 

.016 

.019 

.050 

Normalized  to  Experiment 

3 

1.00 

1.24 

3.20 

5.04 

1.00 

1.19 

3.13 

4.3.4  The  Effects  of  Temperature,  Salinity,  and  Brine  Volume 

Changes  in  temperature  or  salinity  are  physically  manifested  by 
changes  in  brine  volume  according  to  Equation  (3.3).  But,  as  Figure  4.21 
Indicates,  brine  volume  does  not  describe  a  unique  ice  state.  For 
example,  ice  with  a  salinity  of  10  °/0o  and  a  temperature  of  -20°C  has  a 
brine  volume  of  2.8%  as  does  ice  of  salinity  1  °/00  and  temperature  -2°C. 
Their  structure  and  optical  properties,  however,  can  be  quite  different. 

As  discussed  earlier.  Figures  4.17  and  4.18  suggest  that  both  ice 
temperature  and  growth  rate  affect  the  optical  properties  of  sea  ice. 

The  influence  of  temperature  could  best  be  studied  during  the  warming 
phase  of  an  experiment  when  changes  in  ice  thickness,  crystal  structure. 
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and  salinity  were  fairly  small.  Particularly  well  suited  for  this 
purpose,  because  of  the  wide  range  of  temperatures,  was  the  warming 
phase  of  Experiment  6.  Illustrated  in  Figure  4.22  are  the  ice 
temperature  profiles  for  air  temperatures  of  -37,  -22,  and  -2°C. 

The  spectral  albedos  (Figure  4.23)  and  extinction  coefficients 
(Figure  4.24)  are  influenced  quite  strongly  by  ice  temperature. 

Most  dramatic  is  the  drastic  difference  between  ice  above  and  below 
the  eutectic  point.  Below  the  eutectic  the  ice  appears  opaque  and 
white  with  values  of  albedo  and  extinction  coefficient  comparable  to 
those  of  snow.  The  values  sharply  decrease  as  the  ice  temperature 
goes  above  the  eutectic  point.  While  the  change  is  far  less  extreme 
above  the  eutectic  point,  albedo  and  extinction  coefficient  continue 
to  decrease  with  increasing  temperature  as  is  indicated  by  Curves  2 
and  3  in  Figures  4.23  and  4.24.  One  exception  to  this  trend  occurs 
when  water  drains  from  the  ice  surface.  The  surface  texture  then 
becomes  very  crumbly,  resulting  in  an  increase  in  albedo. 

The  data  points  in  Figure  4.25  represent  bulk  albedos  as  a 
function  of  averaged  brine  volume  for  the  warming  phases  of  Experiments 
3,  5,  and  6.  If  all  these  points  are  considered  as  a  whole,  there 
appears  to  be  no  obvious  relationship  between  brine  volume  and  albedo. 
However,  in  each  experiment  the  growth  conditions,  and  hence  the 
distribution  of  brine  within  the  ice,  are  quite  different.  Thus, 
comparisons  should  be  made  only  between  the  individual  experiments. 
Examining  results  from  each  experiment,  a  decrease  in  albedo  with 
increasing  brine  volume  is  evident.  Similarly,  other  results  indicate 
that  the  extinction  coefficient  decreases  with  increasing  brine  volume. 
The  family  of  curves  in  Figure  4.25  also  indicates  that  for  a  given 
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Figure  A. 22;  Temperature  profiles  for  the  wanning  phase  of 
Experiment  6.  The  air  temperature  was:  (1)  -37°C; 

(2)  -22°C;  and  (3)  -2°C. 
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Figure  4.23:  Spectral  albedos  for  the  warming  phase  of  Experiment  6. 

The  numbers  correspond  to  the  same  cases  described  in  Figure  4.22. 

The  cloudy  sky  bulk  albedos  are  (1)  .90,  (2)  .59,  and  (3)  .51. 
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Figure  4.2*:  Spectral  extinction  coefficients  for  the 
warming  phase  of  Experiment  6.  The  numbers  correspond 
to  the  same  cases  described  In  Figure  4.22. 


BRINE  VOLUME 


Figure  4.25:  Brine  volume  versus  bulk  albedo  for  measurements  taken 
during  the  warming  phase  of  (1)  Experiment  6;  (2)  Experiment  5;  and 
(3)  Experiment  3. 
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brine  volume  faster  grown  ice  has  a  larger  albedo. 

The  relationship  of  ice  salinity  to  the  optica]  properties  is  very 
difficult  to  determine,  since  typically  the  effects  of  salinity  are 
intertwined  with  the  effects  of  ice  temperature  and  growth  rate.  The 
influence  of  salinity  is  best  observed  by  comparing  ice  of  different 
salinity,  but  with  similar  thickness,  temperature,  and  growth  history. 
To  investigate  the  role  of  salinity,  two  intermediate  salinity  cases 
(Experiments  7  and  8)  were  performed.  Environmental  conditions  during 
these  two  experiments  were  identical  to  those  of  Experiments  3  and  4, 
except  that  the  water  salinity  was  16  °/00.  With  lower  water  salinity, 
less  brine  was  entrapped  despite  similar  growth  conditions.  By  dupli¬ 
cating  the  air  temperature  of  previous  experiments  comparable  ice 
temperature  profiles  were  obtained.  While  decreased  water  salinity 
did  increase  the  growth  rate  [Equation  (A. 10)]  the  change  was  quite 
small  for  these  experiments. 

Salinity  profiles  and  spectral  albedos  for  sites  of  similar 
thickness  are  plotted  in  Figures  4.26  and  4.27.  Curves  3  and  4  in 
Figure  4.27  comparing  spectral  albedos  from  Experiments  3  and  7,  while 
not  identical,  do  not  display  significant  differences.  However,  there 
is  a  significant  difference  between  Curves  1  and  2,  with  the  less 
saline  ice  of  Experiment  8  displaying  a  greater  albedo.  This  differ¬ 
ence  is  primarily  due  to  a  very  saline  surface  water  skim  present  in 
Experiment  4,  but  not  in  Experiment  8.  The  water  skim  reduces 
scattering  in  the  surface  layer  causing  a  decrease  in  albedo.  Spectral 
extinction  coefficients,  compared  in  Figure  4.28  for  Experiments  4  and 
8,  differ  slightly  for  short  wavelengths  and  are  virtually  identical 
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Figure  4.26:  Salinity  profiles  for  ice  grown  from  water  of  different 
salinities.  In  all  cases  the  ice  is  roughly  22  cm  thick.  The  air 
temperature  and  initial  water  salinity  for  the  cases  are: 

(1)  -20°C,  30.2  °/oo;  (2)  -20°C,  16.5  °/00;  (3)  -10°C,  30.6  °/00; 
and  (4)  -10°C,  16.3  <V00. 
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Figure  A. 27:  Spectral  albedos  for  Ice  grown  from  water  of  different 
salinities.  The  numbers  correspond  to  the  cases  described  in 
Figure  A. 26. 
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WAVELENGTH  (nm) 

Figure  4.28:  Spectral  extinction  coefficients  for  ice 
grown  at  an  air  temperature  of  -20°C  from  water  of 

salinity:  -  30.2  °/00  and - 16.5  °/0o*  The 

numbers  correspond  to  cases  described  in  Figure  4.26. 
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Lane's  (1975)  results  indicate  that  above  5  °/00  salinity  appears 
to  have  little  effect  on  light  transmission.  Certainly  small  amounts 
of  salt  affect  the  optical  properties  of  the  ice;  this  is  quite  evident 
from  the  limiting  case  of  fresh  ice.  However,  for  the  relatively  high 
salinities  investigated  here  and  within  the  resolution  of  this  experi¬ 
ment,  no  consistent  relationship  between  ice  salinity  and  the  optical 
properties  is  evident. 

4.3.5  Discussion 

Since  the  absorption  by  the  ice  and  brine  is  fixed,  changes  in 
the  optical  properties  of  the  ice  result  from  changes  in  the  amount  of 
scattering.  The  initial  sharp  increase  in  albedo  with  ice  thickness 
is  due  to  a  rapidly  increasing  nuiuoer  of  opportunities  for  back  scatter¬ 
ing.  But  as  the  ice  becomes  thicker  it  is  increasingly  difficult  for 
light  to  penetrate  deep  into  the  ice,  be  backscattered ,  and  travel 
back  to  the  surface.  Contributions  to  albedo  are  mostly  from  the 
upper  layers,  hence  the  asymptotic  behavior  with  thickness  as 
illustrated  in  Figure  4.13. 

The  optical  properites,  albedo  in  particular,  are  sensitive  to 
ice  surface  conditions.  When  an  ice  surface  becomes  crumbly  and 
irregular,  the  number  of  air-ice  interfaces  increases,  enhancing 
scattering  and  resulting  in  larger  albedos.  A  water  skim  smooths 
surface  irregularities  reducing  scattering  and  is  a  layer  where 
absorption,  but  no  scattering  takes  place,  effecting  a  decrease  in 
albedo.  This  decrease  is  more  pronounced  at  wavelengths  where  the 
absorption  is  great.  Theoretical  studies  of  how  the  surface  layer 
affects  the  optical  properties  are  discussed  in  Section  5.4.3. 


Growth  rate  affects  the  optical  properties  by  altering  the  ice 
structure.  As  the  growth  rate  increases,  more  brine  is  entrapped 
and  the  platelet  and  crystal  sizes  decrease  (Weeks  and  Hamilton,  1962; 
Weeks  and  Assur,  1967).  Weeks  and  Assur  (1967)  empirically  related 
platelet  size  and  growth  rate  by 

f  =  -2-  (4.12) 

a2 

where  a  is  the  platelet  width  (cm),  f  is  the  growth  rate  (cm/sec),  and 
p  is  1.33  x  10-7(cm3/sec) .  Predicted  platelet  widths,  at  a  depth  of 
10  cm,  range  from  .28  mm  for  Experiment  6  to  .54  mm  for  Experiment  3. 
Smaller  platelets  and  crystals  imply  more  scattering  locations  per  unit 
volume.  Crystal  size  and  platelet  width  also  gradually  increase  with 
depth  (Weeks  and  Hamilton,  1962) .  Crystals  and  platelets  in  the  tank 
experiments,  while  not  rigorously  monitored,  appeared  to  follow  these 
trends.  Examining  air  bubbles  in  fresh  ice,  Gow  and  Langston  (1972) 
found  as  growth  rate  increases  the  average  bubble  radius  decreases  and 
the  density  of  bubbles  increases.  Higher  growth  rates  yielded  more 
brine  inclusions,  platelets,  crystals,  and  air  bubbles,  producing  more 
scattering  sites.  This  increase  in  the  number  of  scattering  sites 
produces  increases  in  albedo  and  extinction  coefficient. 

Below  the  eutectic  point  ice  appearance  dramatically  changes, 
becoming  very  white  with  solid  salts  visible.  Sinha  (1977),  using  an 
electron  microscope,  photographed  brine  pockets  showing  precipitated 
solid  salts  at  a  temperature  of  -30°C.  The  precipitation  of  solid 
salts,  providing  additional  scattering  locations,  causes  a  marked 
increase  in  albedo  and  extinction  coefficient. 


Above  the  eutectic  point  let*  structure  variations  with  tempera!  me 
arc  less  obvious.  Viewing  a  vertical  lee  sect  ion  in  ttansmiiteii  light, 
as  the  temperature  is  raised,  the  features  within  the  ice  become  more 
visible  and  1 tght  transmission  increases,  but  the  structural  cause  is 
not  apparent  without  detailed  microscopic  examination. 

After  growing  an  Ice  sheet  roughly  20  cm  thick  at  an  ait  tempera 
ture  of  -30°C,  a  core  was  removed  and  a  sample  was  prepared  lor 
observation  under  the  microscope.  The  sample  was  a  horizontal  thin 
section,  roughly  .7  cm  thick  from  10  cm  deep  in  the  Ice.  the  salinity 
of  the  sample  was  11.2  °/OD.  At  this  depth  the  crystals  were  orderly 
and  oriented  with  horizontal  c-axls,  and  a  horizontal  sect  Ion,  cut 
perpendicular  to  the  orientation  ot  the  crystals  and  platelets  includes 
many  boundaries.  The  sample  was  frozen  on  a  glass  plate  and  lurthet 
thinned  using  a  bandsaw;  finally  its  surface  was  polished.  To  minimize 
brine  drainage  from  the  sample  the  preparation  was  performed  at  a 
temperature  of  -25°C.  Finally  water  was  frozen  around  the  periphery 
of  the  sample  to  seal  the  edges.  Inhibiting  brine  drainage.  the 
sample  was  placed  under  the  microscope  and  warmed  over  a  period  ot  two 
days.  The  procedure  was  to  increase  the  air  temperature,  allow  the 
ice  to  attain  thermal  equilibrium,  photograph  the  sample,  and  repeat 
the  process.  This  was  done  at  fourteen  temperatures  between  2  7‘'t'  and 
-1°C. 

Four  microphotographs,  taken  at  temperatures  of  -7,  4,  and 

-1°C  are  presented  In  Figure  4,29  a-d.  Arrows  are  vised  to  follow  the 
evolution  of  representative  features  as  the  brine  volume  increases 
with  temperature.  The  thin  dark  lines  between  platelets  in  Figure  4..’*>A 


Figure  4.29:  Microphotographs  of  a  horizontal  thin  section:  (A) 
(D)  — 1 0 C .  The  double  arrow  in  Figure  4.29A  is  1  mm  in  length. 
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are  due  to  the  precipitated  salts.  These  lines  do  not  appear  in  the 
photographs  (B,  C,  and  D)  taken  above  the  eutectic  point. 

The  microscope  observations  illustrate  several  points  regarding 
changes  in  ice  structure  during  warming.  As  the  brine  volume  increases, 
the  platelet  boundaries  become  broader  and  smoother.  Small  isolated 
pockets  of  brine  enlarge  and  sometimes  connect.  Similarly,  small  air 
bubbles,  released  within  the  ice  by  melting  often  join  together. 

These  changes  all  tend  to  reduce  the  amount  of  scattering  within  the 
ice,  and  thus  reduce  the  albedo  and  the  extinction  coefficient. 

To  avoid  confusion  when  examining  the  effects  of  brine  volume  on 
the  optical  properties,  the  influences  of  warming  and  growth  rate  must 
be  distinguished.  For  two  ice  sheets,  grown  at  different  rates  and 
brought  to  the  same  temperature,  the  faster  grown  saltier  ice  will 
have  a  larger  brine  volume;  but,  since  it  has  a  greater  number  of  brine 
inclusions,  scattering  is  greater  and  the  albedo  and  extinction 
coefficient  larger.  Dichotomies  regarding  the  relationship  between 
brine  volume  and  the  optical  properties  are  resolved  by  careful 
consideration  not  only  of  the  brine  volume,  but  of  the  brine 


distribution. 


CHAPTER  FIVE 


FOUR  STREAM  MODELING 


5.1  INTRODUCTION 

To  this  point  the  optical  properties  of  sea  ice  have  been  examined 
using  a  Dunkle  and  Bevans  photometric  model.  To  introduce  realistic 
scattering  functions  and  investigate  the  physics  of  the  radiative 
transfer  the  four  stream  DOM  model  was  used  to  examine  the  influence  of 
,  optical  depth,  and  phase  function  on  albedo  and  transmittance. 
Experimental  results  were  fitted  by  one  or  two  layer  models.  Using 
suitable  values  for  the  optical  parameters  the  effects  of  wavelength, 
ice  thickness,  incident  radiation  field,  and  refraction  were  determined. 

Calculations,  for  the  most  part,  were  performed  at  650  nm.  This 
particular  wavelength  was  selected  since  the  spectrophotometer  had  a 
large  signal  to  noise  ratio  and  reasonably  high  resolution.  In  addition, 
it  was  possible  to  use  phase  functions  obtained  from  preliminary  results 
of  scattering  experiments  conducted  with  a  HeNe  laser  (Grenfell  and 
Hedrick,  personal  communication,  1978).  In  a  few  cases  the  results  at 
650  nm  were  extended  to  calculate  spectral  albedos  and  transmi ttances 
from  400  to  1000  nm.  Except  for  the  direct  beam  cases  discussed  in 
Section  5.4,  an  isotropic  incident  radiation  field  was  used  for  all 


calculations . 


5.2  PARAMETER  STUDY 


5.2.1  Parameters  of  the  Discrete  Ordinates  Method 

The  parameters  of  the  discrete  ordinates  method  are  five  quantities: 
the  single  scattering  albedo  (wq),  the  optical  depth  (i^),  the  scattering 
coefficient  (0),  the  absorption  coefficient  (k) ,  and  the  depth  (H) ,  which 
are  related  by  Equations  (2.18)  and  (2.22).  The  phase  function  [p(0)] 
and  the  incident  radiation  field  must  also  be  specified  to  calculate 
radiances  and  irradiances  as  a  function  of  depth. 

Although  k,  o,  and  H  are  the  fundamental  physical  properties  of  the 
medium,  it  is  often  more  convenient  to  discuss  dip  and  T^,  which  are  the 
explicit  parameters  appearing  in  the  theory.  In  practice,  the  absorption 


This  relationship  holds  true  at  any  depth  within  the  ice.  For  a 
particular  value  of  uig  physical  and  optical  thickness  are  linearly 
related,  with  the  slope  increasing  for  larger  uiQ  (Figure  5.1a).  For  the 
four  values  of  illustrated,  .97,  .98,  .99,  and  .995,  an  optical 
thickness  of  one  corresponds  to  a  physical  thickness  of  9.38,  6.25,  3.13, 
and  1.56  cm  respectively.  Figure  5.1b  shows  the  steepening  increase  in 
with  increasing  iLq  at  a  given  physical  thickness.  Figure  5.1  is 
useful  in  relating  the  perceptual  physical  thickness  of  the  ice  to  the 
conceptual  optical  thickness  of  the  model. 

Sea  ice,  being  a  medium  dominated  by  scattering,  has  large  values 
of  uig.  For  the  salt  ice  studied  in  the  tank  experiments,  u'0  typically 
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ICE  THICKNESS  (cm) 


Figure  5.1:  Interrelation  between  single  scattering  albedo  (d tp), 
optical  thickness  (x^) ,  and  physical  (ice)  thickness  (H) 

(k  ■  .0032  cm-1).  (A)  Physical  thickness  versus  optical 

thickness  for  constant  values  of  w0:  (1)  d>0  -  .97;  (2)  u>p  -  .98; 

(3)  dip  >»  .99;  and  (4)  u>0  ■  .995.  (B)  Optical  thickness  versus 

w0  for  constant  values  of  H:  (1)  H  ■  5  cm;  (2)  H  ■  15  cm;  and 
(3)  H  -  25  cm. 
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ranged  from  .95  to  .995,  H  from  0  to  30  cm,  and  i„  from  0  to  19.  While 
determines  the  amount  of  scattering,  it  is  the  phase  function  which 
specifies  the  angular  dependence  of  scattering.  Commonly  used  in  studies 
of  radiative  transfer  is  the  Henyey-Greenstein  phase  function 

1  _  o- 

p(0)  -  - - - 8 - —  (5.2) 

(l+g"-2g  cos 

where  the  asymmetry  factor,  g,  determines  the  relative  importance  of 
forward  and  backward  scattering.  The  scattering  is  isotropic  for  g  =  0 
and  becomes  more  forward  peaked  as  g  approaches  one.  Mathematically 
Henyey-Greenstein  phase  functions  are  convenient  since  they  yield  a 
simple  relationship  for  the  Legendre  expansion  coefficients 

-  (2m)gc  .  (5.3) 

Examples  of  the  Henyey-Greenstein  phase  function  for  g  equal  to  0,  .5, 
and  .75  are  plotted  in  Figure  5.2. 

Grenfell  and  Hedrick  (personal  communication,  1978),  using  a  HeNe 
laser  as  a  light  source,  measured  the  angular  dependence  of  scattering 
by  ice  samples  from  Experiments  3  and  A.  Figure  5.3,  Curve  1  is  an 
eyeballed  best  fit  of  their  preliminary  results.  The  importance  of 
backscattering  is  indicated  by  the  increase  of  the  phase  function  for 
angles  larger  than  105°.  Because  of  this  increase,  a  one  parameter 
Henyey-Greenstein  phase  function  was  insufficient  to  describe  the 
scattering.  Used  instead  was  a  three  parameter  Henyey-Greenstein  phase 


function  of  the  form 


PHASE  FUNCTION 


Figure  5.3:  Phase  function  determination  for  scattering  experiment 
results  (Grenfell  and  Hedrick,  personal  communication): 

(1)  experimental  results;  (2)  best  fit  using  a  three  parameter 
Henyey-Greenstein  phase  function;  and  (3)  phase  function 
calculated  using  w^'s  and  Equation  2.21. 
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p(0)  =  (1-y) 


1_8l 


1- 


+  y- 


®2 


(l+gj-ZgjCos  6)3/2  (l+g2-2g2cos  0 )  3  ^  2 


(5. A) 


where  the  second  term  enhances  the  backscatter ing  (g2  <  0)  and  y  is  the 
weighting  factor.  A  nonlinear  regression  routine  using  data  digitized 
at  fifteen  degree  intervals  was  used  to  determine  the  three  c  fflcients. 
Best  fit  coefficients  of  y  =  .279,  gi  =  .71,  and  g2  =  -.25  were 
substituted  into  Equation  (5.4),  which  was  plotted  as  Curve  2. 

Introducing  this  phase  function  into  Equation  (2.23),  the  associated 
w^'s  were  determined  by  numerical  integration,  then  substituted  into 
Equation  (2.21)  to  obtain  Curve  3.  The  four  term  expansion  of  this  phase 
function  is  quite  good  at  *U2  and,  though  less  precise  at  is 

accurate  for  calculations  of  irradiance. 


5.2.2  Sensitivity  of  Albedo  and  Transmittance  to  the  DOM  Parameters 
The  Phase  Function 

Using  the  single  parameter  Henyey-Greenstein  phase  function 
described  by  Equation  (5.2),  albedos  and  transmittances  were  calculated 
over  a  range  of  g's.  Values  at  650  nm  for  ice  20  cm  thick  and  various 
(lip's  are  plotted  in  Figure  5.4.  The  largest  albedos  and  smallest 
transmittances  are  produced  by  isotropic  scattering  (g  =  0) .  As  the 
asymmetry  factor  (g)  increases,  and  the  scattering  becomes  more  forward 
peaked,  the  albedo  decreases  and  the  transmittance  increases. 

Thickness  and  Single  Scattering  Albedo 

Isopleths  of  albedo  and  transmittance  for  pertinent  ranges  of 
and  H  are  plotted  in  Figures  5.5  and  5.6.  The  calculations  were 


ALBEDO 
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ASYMMETRY  FACTOR 


Figure  5.4:  Effect  of  asymmetry  parameter  g  on  albedo  (— )  and 

transmittance  (■ - ).  Calculations  are  for  ice  20  cm  thick 

at  a  wavelength  of  650  nm  using  the  three  parameter  Henyey- 
Greenstein  phase  function  with  diffuse  incident  radiation: 
(1)  u>0  =  .96;  (2)  w0  =»  .98;  and  (3)  =  .99. 


TRANSMITTANCE 


Figure  5.5:  Isopleths  of  albedo  as  a  function  ol 
thickness,  calculated  at  650  nm  using  a  three 
Henyey-Greenstein  phase  function  with  diffuse 
radiat ion . 
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performed  at  650  nm  using  the  three  parameter  Henvey-Greenstein  phase 
function.  For  a  given  dig ,  as  the  Ice  thickness  increases,  the  albedo 
increases  and  the  transmittance  decreases.  Similarity,  holding  ice 
thickness  constant  while  increasing  effects  an  increase  in  albedo 
and  a  decrease  in  transmittance.  As  u>q  increases  albedo  and  trans¬ 
mittance  become  more  sensitive  to  changes  in  either  ,1^  or  H. 

The  Effects  of  Refraction 

Using  an  appropriate  phase  function  [Equation  (5.4)1  and  reasonable 
values  for  and  depth,  albedos  and  transmittances  for  both  nonrefract¬ 
ing  and  refracting  media  were  calculated  under  direct  beam  and  diffuse 
incident  conditions.  The  angle  of  the  direct  beam  was  representative 
for  arctic  summer  conditions.  As  illustrated  in  Table  5.1,  in  all  cases, 
refraction  caused  a  decrease  in  albedo  and  an  increase  in  transmittance. 

In  the  case  of  refraction  a  portion  of  the  upwelling  light  at  the 
interface  is  totally  internally  reflected  downwards  into  the  medium, 
thus  decreasing  the  albedo  and  increasing  the  transmittance.  Kefraction 
converges  the  incident  radiation  into  a  more  normal  direct  ion  also 
increasing  transmission. 

5.3  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  ALBEDOS  AND  TRANSMITTANCES 

AT  650  nm 

5.3.1  One  Layer  Model 

The  four  stream  model  was  tested  using  results  from  the  tank 
experiments.  While  ice  thickness  data  were  available  for  each  experi¬ 
mental  site,  single  scattering  albedos  and  phase  functions  are  unknown. 


The  three  parameter  Henyey-Greensteln  phase  function  I  Equation  (S.41I 
discussed  In  Section  5.2.1  was  assumed  to  describe  t lie  angular 
scattering  for  all  cases.  Values  of  u>(1  were  determined  using  an 
inverse  procedure.  For  a  particular  experimental  site  was  selected 
so  that  the  albedo  calculated  using  the  four  stream  model  was  equal  to 
the  experimental  albedo.  The  analysis  was  done  at  a  wavelength  of 
b50  nm.  In  order  to  compensate  for  the  effects  of  instrumental  band¬ 
width  on  the  results,  an  absorption  coefficient  oi  .0042  cm'1  was  used 
rather  than  the  reported  value  at  650  nm  of  .0032  cm"'  (Sauberer,  1950). 
This  value  of  .0042  cm"'  was  determined  by  convolving  a  tvpical  spectral 
transmission  curve  with  the  bandpass  of  the  instrument . 

This  technique  applied  to  data  from  eleven  experimental  sites  and 
the  results  summarized  in  Table  5.2.  Theoretical  and  experimental 
albedos  were  matched  quite  well  by  values  of  ranging  from  .955  to 
.9997.  However,  agreement  between  theoretical  and  experimental  trans¬ 
mit  tances  is  not  as  good.  In  all  cases  except  two,  the  theoretical 
transmittance  is  larger  than  the  experimental,  with  differences  ranging 
from  6  to  55X. 

There  are  a  number  of  possible  explanations  for  the  lack  of  agree¬ 
ment  between  theoretical  and  experimental  transmi t tances .  The  assump¬ 
tion  of  homogeneous  ice,  implicit  in  a  one  layer  model,  is  probably  the 
source  of  error  in  the  theoretical  calculations.  Temperature  gradients, 
changes  in  growth  rates,  and  differences  in  crystal  orientation  can  all 
cause  vertical  variations  in  the  physical  state  of  the  ice,  which  affect 
the  optical  properties.  In  addition,  changes  in  the  ice  surface  strongly 
influence  the  albedo  and  transmittance.  Table  5.2  indicates  that  the 
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closest  agreement  between  theoretical  and  experimental  transraittances 
occurs  when  the  ice  is  most  nearly  homogeneous.  The  difference  is 
smallest  for  thin  ice  and  warm  ice.  Conversely,  the  largest  differences 
are  exhibited  for  the  very  inhomogeneous,  cold,  fast  grown  ice.  A  one 
layer  model  thus  appears  to  be  inadequate  for  describing  the  transfer 
of  radiation  in  vertically  complex  ice. 

5.3.2  Two  Layer  Model 

In  an  effort  to  obtain  better  agreement  between  theoretical  and 
experimental  results,  a  two  layer  model  was  formulated.  Aside  from 
the  boundary  conditions  the  theoretical  development  of  a  two  layer 
model  was  the  same  as  the  one  layer  model.  To  calculate  the  constants 
of  integration,  Lj's,  in  a  two  layer  model,  a  system  of  eight  equations 
and  eight  unknowns  must  be  solved.  The  four  additional  boundary 
conditions  are  determined  by  requiring  continuity  of  radiance  at  the 
interface  between  the  two  layers 


X1 (~yl »Tl) 

=  I2(-y1,x1) 

(5.5) 

1 1  (“1*2  ♦ T  1  > 

(5.6) 

I] (M  j »Tj) 

=  ^2  ^1  *T  1  ^ 

(5.7) 

1 1  T  1 ) 

=  ^2  ^2  *  ^  1 ) 

(5.8) 

where  tj  is  the  optical  depth  of  layer  one  and  Ij  and  I2  are  the 
radiances  in  layers  one  and  two. 

The  two  layer  model  was  applied  to  the  results  from  Experiments  5 
and  6.  The  two  cases  selected  had  the  largest  temperature  gradients 
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and  were  poorly  represented  by  the  one  layer  model.  This  is  particularly 
true  of  Experiment  6,  where  a  few  centimeter  thick  surface  layer  was 
below  the  eutectic  point,  resulting  in  a  region  containing  precipitated 
salts  which  was  very  different  in  appearance  from  the  underlying  ice. 

Using  realistic  values  of  <I>0  and  layer  thickness  a  range  of  two 
layer  cases  were  examined.  The  two  which  best  represent  the  experimental 
results  are  summarized  in  Table  5.3.  While  the  agreement  between 
experimental  and  theoretical  results  is  closer  for  a  two  layer  model, 
substantial  differences  still  remain. 


Table  5.3 


TWo  Layer.  Four  Stream  Modeling  of  Cold  Ice 


LAYER  I 

LAYER  2 

THEORETICAL 

EXPERIMENTAL  | 

Experiment 

Air  Temp 
°C 

.  Thickness 

“o  (cm) 

, 

thickness 

(c») 

a 

T* 

a 

T* 

6 

-37 

.99993  2.0 

.97 

24.0 

.898 

.018 

.959 

.002 

5 

-30 

.992  11.0 

.98 

13.0 

.599 

.282 

.587 

.147 

"Transmitted  at  19  cm. 


A  more  detailed  four  stream  analysis  of  the  experimental  results 
is  not  warranted  until  better  data  concerning  values  of  the  single 
scattering  albedo  and  the  phase  function  become  available.  Experiments 
explicitly  measuring  w0  and  the  phase  function  for  a  variety  of  ice 
types  and  conditions  are  necessary.  Only  then  can  multilayered  models, 
accurately  representing  ice  surface  and  interior  conditions,  be 
implemented. 
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5.4  FURTHER  APPLICATIONS  OF  THE  FOUR  STREAM  MODEL 

Modeling  of  the  experimental  results  indicates  that  for  the  three 
parameter  Hcnvev  Greenstein  phase  function  and  l Equation  (5.41],  uy, 
ranges  from  .*155  for  warm  slowly  grown  ice  to  as  high  as  .99993  for  ice 
below  the  eutectic  point.  Using  appropriate  parameter  values  the 
four  stream  model  is  now  extended  to  theoretically  investigate 
additional  topics. 


5.4.1  Depth  Dependence  of  Irradiance 

Though  the  focus  of  the  four  stream  modeling  has  been  on  the 
calculation  of  albedos  and  transmi ttances ,  depth  dependent  upwelling, 
downwelling,  and  net  irradiances  can  also  be  determined.  Figures  5.7 
and  5.8  illustrate  normalized  downwelling  and  net  irradiance  as  a 
function  of  depth  at  650  nm.  The  irradiances  were  calculated  for  a 
single  layered  ice  slab,  25  cm  thick,  under  diffuse  illumination  using 
a  three  parameter  Henvey-Greenstetn  phase  function  and  values  of  single 
scattering  albedo  of  .97,  .98,  and  .99.  Evident  from  Figure  5.7  is  a 
decrease  in  downwelling  irradiance  with  depth,  a  decrease  which  becomes 
more  pronounced  as  increases.  In  addition,  for  a  given  ioq  the  slope 
of  F+(z)  becomes  steeper  with  depth  (d'F4 (z) /dz-  < 0) .  The  slope  of  the 
net  irradiance,  the  Important  quantity  when  considering  the  energy 
absorbed  in  the  ice,  decreases  as  depth  increases.  In  general,  the 
extinction  coefficient  is  expressed  as 


F  dz  dz 


(5.9) 


DEPTH  (cm) 


Figure  5.7:  Normalized  downwelllng  lrradiance  at  650  nm  as  a  function 


of  depth  for:  (1)  J>0  =  .97;  (2)  -  . 

Calculations  were  performed  for  diffuse 


98;  and  (3)  ”  .99. 

Incident  radiation  using 


the  three  parameter  Henyey-Greenstein  phase  function. 


DEPTH  (cm) 


Figure  5.8:  Normalized  net  irradiance  at  t>50  nm  as  a  function  of  depth 
for:  (1)  u>0  »  .97;  (2)  d»0  -  .98;  and  (3)  ^  -  .99.  The  three 

parameter  Henyev-Greensteln  phase  function  was  used  and  the  incident 
radiation  field  was  diffuse. 
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where  Is  the  net  irradianee  at  depth  z  (Grenfell  and  Maykut ,  1978). 

Thus  the  slope  of  the  curves  in  Figure  5.8  is  the  extinction  coefficient, 
which  asymptotically  approaches  a  minimum  value  as  the  depth  increases. 
The  changing  slope  demonstrates  the  inapplicability  of  the  Bouguer- 
Lambert  law  in  this  region. 

5.4.2  Direct  Beam  Incident  Radiation 

When  using  the  discrete  ordinates  method  it  is  not  necessary  to 
assume  isotropic  Incident  radiation.  Clear  sky  conditions  can  also  he 
investigated  by  representing  the  incident  radiation  field  by  a  direct 
beam  at  -p . 

Considering  a  one  layer  ice  slab  25  cm  thick,  the  effect  of  the 
angle  of  incidence  on  albedo  and  transmittance  for  clear  sky  conditions 
is  displayed  in  Figure  5.9.  Three  values  of  ic.  were  examined  with  a 
three  parameter  Henyey-Greenstein  phase  function  used  in  each  case. 

Beyond  80°  albedo  and  transmittance  are  very  sensitive  to  changes  in  the 
angle  of  incidence.  This  is  a  result  of  the  increase  in  the  Fresnel 
reflection  of  the  direct  beam.  At  approximately  54°,  results  for  direct 
beam  incident  radiation  are  equal  to  those  for  diffuse  incident  radiation. 


ALBEDO 
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Figure  5.9:  Albedo  ( - )  and  transmittance  ( - )  as  a  function 

of  solar  angle  (angle  of  direct  incident  beam)  for:  (1)  iop  *  .97; 
(2)  w  -  .98;  and  (3)  do0  -  .99  (values  for  isotropic  incident). 

The  three  parameter  Henyey-Greenstein  phase  function  was  used  in 
cloudy  sky  results. 
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5.4.3  Spectral  Albedos  and  Transmittances 

By  using  a  modified  version  of  the  four  stream  model  spectral 
albedos  and  transmittances  were  calculated  at  50  nm  intervals  between 
400  and  1000  nm.  To  accomplish  this  it  was  necessary  to  determine  the 
wavelength  dependence  of  the  basic  four  stream  parameters:  the  absorp¬ 
tion  coefficient,  the  scattering  coefficient,  and  the  phase  function. 
Summarized  in  Table  5.4  are  spectral  absorption  coefficients  measured 
by  Sauberer  (1950)  for  pure  bubble  free  ice.  With  scattering  inhomo¬ 
geneities  much  larger  than  optical  wavelengths  geometrical  scattering 
dominates  in  sea  ice.  Since  the  index  of  refraction  is  a  weak  function 
of  wavelength,  the  scattering  coefficient  and  the  phase  function  are 
assumed  to  be  independent  of  wavelength  from  400  to  1000  nm.  There  is, 
however,  no  experimental  verification  of  this  independence,  though 
investigations  are  presently  being  conducted  (Grenfell,  personal 
communication) . 

Theoretically  calculated  spectral  albedos  are  compared  to  experi¬ 
mental  values  in  Figure  5.10.  Ice  21.1  cm  thick  grown  at  an  air 
temperature  of  -10°C,  was  selected  for  comparison  since  the  ice  is 
fairly  homogeneous  and  thus  suitable  for  one  layer  modeling.  Using  a 
value  for  ui()  of  .974  (Table  5.2),  a  scattering  coefficient  of  0.119  is 
calculated.  Theoretical  and  experimental  albedos  closely  agree  below 
800  nm,  but  above  this  experimental  values  are  somewhat  higher.  The 
disagreement  in  this  region  is  probably  due  to  the  wide  bandpass  of  the 
instrument . 

The  magnitude  and  spectral  dependence  of  albedo  and  transmittance 
are  markedly  influenced  by  ice  thickness.  Using  a  scattering  coefficient 
of  .157  cm-1,  which  at  650  nm  corresponds  to  a  utn  of  .98,  spectral 


Table  5.4 

Absorption  Coefficients  of  Ice 
(Snubeier ,  10^0) 


X  (nra) 

k(cm- 1 ) 

400 

.0003 

450 

.0005 

500 

.0008 

550 

.0013 

6C0 

.0020 

650 

.0032 

700 

.0060 

750 

.0110 

800 

.0200 

850 

.0330 

900 

.0590 

950 

.1050 

1000 

.2500 

ALBEDO 
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Figure  3.10:  Comparison  of  spectral  albedos  for  lee  21  cm  thick  at 

an  air  temperature  of  -10°C:  ( - )  experimental  and  ( - ) 

theoretical.  Theoretical  calculations  used  the  three  parameter 
Henyey-Greensteln  phase  function  and  diffuse  incident  conditions. 
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albedos  (Figure  5.11)  and  transmit tnnces  (Figure  5.12)  were  calculated 
for  Ice  ranging  in  thickness  from  5  to  25  cm.  As  the  ice  thickens,  the 
increase  in  albedo  is  greater  at  shorter  wavelengths  where  u>0  is  larger. 
This  results  in  a  wavelength  dependence  for  albedo  which  is  weak  for 
thin  ice  and  becomes  more  pronounced  as  the  thickness  increases.  The 
transmittance  sharply  decreases  at  longer  wavelengths  and  as  the  ice 
thickness  increases,  narrowing  the  spectrum  of  transmitted  light. 

The  spectral  albedo  of  an  ice  cover  is  affected  by  the  nature  of 
the  surface  layer.  To  investigate  the  extent  of  this  effect  spectral 
albedos  for  a  two  layer  model  are  calculated.  Five  different  cases 
were  studied  (Figure  5.13).  Curve  3  represents  25  cm  thick  ice  with  a 
scattering  coefficient  of  .157.  In  Curves  1  and  2  the  top  five  and  one 
centimeters  of  ice  were  replaced  by  water  of  negligible  scattering 
coefficient.  It  was  assumed  that  the  water  layer  did  not  change  the 
scattering  properties  of  the  underlying  ice.  This  is  somewhat  idealized 
since  the  scattering  due  to  the  surface  irregularities  of  the  ice  Is 
reduced  by  the  presence  of  water.  Curves  4  and  5  demonstrate  the  effect 
of  highly  scattering  (o  «  3.997  cm-1),  one  and  five  centimeter  thick 
surface  layers,  representative  of  ice  below  the  eutectic  point. 

The  presence  of  a  surface  water  layer  decreases  the  albedo  while 
the  albedo  is  enhanced  by  a  surface  scattering  layer.  The  thicker  the 
surface  layer  the  more  pronounced  the  effect.  The  influence  of  a  surface 
scattering  layer  is  dependent  on  wavelength,  as  is  displayed  by  the 
change  in  shape  of  the  spectral  albedo  curve.  As  has  been  indicated 
earlier  (Figures  5.5  and  5.11),  the  thickness  of  ice  significantly 
contributing  to  the  albedo  decreases  with  d>.  (for  the  limiting  case  of 
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Figure  5.11:  Spectral  albedos  calculated  for  Ice  ranging  In  thickness 
front  2.5  to  25  cm  using  a  scattering  coefficient  of  .157.  The 
three  parameter  Henyey-Greenstein  phase  function  and  diffuse 
incident  conditions  were  used. 


TRANSMITTANCE 


128 


400  600  800  1000 

WAVELENGTH  (nm) 


Figure  5.12:  Spectral  transmit tances  calculated  for  ice  ranging  in 
thickness  from  2.5  to  25  cm  using  a  scattering  coefficient  of  .157. 
The  three  parameter  Henyey-Greenstein  phase  function  and  diffuse 
incident  radiation  were  used. 
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WAVELENGTH  (nm) 


Figure  5.13:  Effect  of  surface  layer  on  albedo.  The  total 
thickness  is  25  cm  and  the  scattering  coefficient  for  the 
underlying  ice  is  .157:  (1)  5  cm  thick  surface  water  layer; 

(2)  1  cm  thick  surface  water  layer;  (3)  no  surface  layer; 

(4)  1  cm  thick  surface  scattering  layer  (o  ■  3.997);  and 

(5)  5  cm  thick  surface  scattering  layer  (o  *=  3.997).  Diffuse 

incident  conditions  and  the  three  parameter  Henyey-Creenstein 
phase  function  were  used. 
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£>0  =  0  the  albedo  is  due  only  to  specular  reflection)  .  Therefore  the 
albedo  at  longer  wavelengths,  with  larger  absorption  coefficients,  is 
more  sensitive  to  surface  conditions. 


CHAPTER  SIX 


SUMMARY  AND  SUGGESTIONS  FOR  FUTURE  RESEARCH 

The  optical  properties  of  cold,  young  ice  are  markedly  different 
from  those  of  melting  first-year  ice  and  older  multiyear  ice.  In 
particular,  extinction  coefficients  can  be  as  much  as  15  times  larger. 
Under  cold  conditions  albedos  for  thin  ice  (25  cm)  are  comparable  to 
values  for  older,  thicker  ice. 

While  these  experiments  do  not  allow  a  quantitative  statement  of 
relationships  between  the  optical  properties  and  the  physical  parameters 
of  the  ice,  some  generalizations  can  be  made.  For  young  growing  ice  the 
albedo  increases  as  the  ice  thickens.  The  albedo  is  quite  sensitive  to 
ice  surface  conditions,  with  a  crumbly  surface  causing  an  increase  in 
albedo,  a  water  skim  a  decrease,  and  the  presence  of  solid  salts  a 
dramatic  increase  to  values  comparable  to  snow.  When  comparing  ice  of 
similar  temperature  profiles  and  thicknesses,  but  with  different  growth 
histories,  the  faster  grown  ice  has  larger  albedos  and  extinction  coeffi¬ 
cients.  Over  the  range  of  salinities  investigated  (4  °/OD  to  14  °/00) 
ice  salinity  does  not  discernibly  influence  the  optical  properties.  As 
the  temperature  of  the  ice  decreases,  its  albedo  and  extinction  coeffi¬ 
cient  increase.  When  the  ice  temperature  drops  below  the  eutectic  point 
the  increase  is  discontinuous  and  dramatic  due  to  the  presence  of 
precipitated  salts.  The  brine  volume  of  the  ice  changes  with  temperature 
and  thus  is  related  to  the  optical  properties.  When  comparing  brine 
volumes  of  ice  grown  under  different  conditions,  the  influence  of  brine 
volume  is  less  clear  since  other  parameters  such  as  growth  rate  are  also 
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significant.  It  is  important  to  consider  the  distribution  of  brine  as 
well  as  the  brine  volume  when  relating  it  to  the  optical  properties. 

Changes  in  ice  structure  which  enhance  scattering  result  in 
increased  albedos  and  extinction  coefficients.  A  crumbly  surface  with 
many  air-ice  interfaces  is  highly  scattering,  while  a  water  skim  with 
water-ice  interfaces  is  not.  With  an  index  of  refraction  of  1.52  and 
small  dendritic  and  cubic  crystals,  solid  salts  deposited  on  the  ice 
surface  optically  act  similar  to  a  thin  snow  layer.  Faster  grown  ice 
has  a  complex  structure  with  smaller  and  more  disordered  platelets  and 
crystals.  As  the  ice  warms  and  the  brine  volume  increases,  edges  become 
smoother,  brine  pockets  expand,  and  air  bubbles  coalesce,  all  of  which 
decrease  scattering. 

Future  laboratory  work  should  focus  on  this  problem  of  relating  the 
optical  properties  to  the  physical  state  and  structure  of  the  ice.  Much 
of  the  previous  work  regarding  the  optical  properties  of  sea  ice  has 
been  primarily  phenomenological.  Albedos  and  extinction  coefficients 
were  determined  and  associated  with  various  categories  of  ice  while  the 
radiative  processes  occurring  in  the  ice  were  largely  ignored.  In  this 
research,  by  combining  observations  of  changes  in  the  optical  properties 
as  the  physical  state  of  the  ice  varied  with  a  more  advanced  photometric 
model,  some  qualitative  insight  into  radiative  transfer  in  sea  ice  was 
obtained. 

Complex  enough  to  accurately  describe  the  physics  involved,  yet 
simple  enough  to  avoid  time  consuming  computing  techniques,  the  four 
stream  model,  considering  both  anisotropic  scattering  and  refraction  at 
the  boundaries,  provides  an  effective  method  to  examine  the  radiative 
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processes  In  sea  ice.  Present  modeling  efforts  are  severely  limited  by 
a  scarcity  of  data  concerning  the  scattering  coefficients  and  the  phase 
functions . 

Most  productive  would  be  a  series  of  experiments  on  small  homogeneous 
ice  samples  to  determine  single  scattering  albedos  and  phase  functions 
over  an  extensive  range  of  conditions.  Unlike  the  large  heterogeneous 
ice  slabs  of  the  tank  experiments,  these  small  samples  could  be  physi¬ 
cally  uniform,  allowing  a  more  precise  statement  of  the  relationship 
between  the  physical  parameters  and  the  single  scattering  albedo  and 
phase  function.  The  scattering  properties  of  a  given  sample  could  he 
measured  as  a  function  of  temperature,  wavelength,  and  sample  orientation. 
Measurements  of  different  samples  could  indicate  the  effects  of  salinity, 
platelet  size,  growth  conditions,  air  bubble  density  and  size  distribu¬ 
tion,  and  depth  within  the  ice.  In  this  fashion  the  optical  influence 
of  specific  physical  parameters  could  be  isolated.  Also  necessary  are 
more  accurate  measurements  of  the  spectral  absorption  coefficients  of 
pure  bubble  free  ice,  particularly  at  short  wavelengths  where  the 
absorption  is  very  small. 

Calculations,  using  eight  and  sixteen  stream  discrete  ordinate 
models,  are  needed  to  investigate  the  angular  distribution  of  upwelling 
radiance.  At  present,  when  considering  refraction  at  the  ice-air 
interface,  the  upwelling  radiance  is  defined  to  be  a  linear  function  of 
p.  Examining  a  few  representative  cases,  the  accuracy  of  this  assumption 
could  be  further  tested. 

This  additional  information  would  enable  a  more  accurate  and 
extensive  photometric  modeling  of  sea  ice.  Scattering  results  from 
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small  homogeneous  samples  could  be  assembled  using  a  multilayer  model 
to  optically  represent  a  very  complex  ice  cover.  Changes  in  the  deposi¬ 
tion  of  radiant  energy  within  an  evolving  ice  cover  could  be  examined. 

By  coupling  experimentally  determined  phase  functions  and  single 
scattering  albedos  with  four  stream  photometric  models,  an  enhanced 
understanding  of  radiative  transfer  processes  in  sea  ice  could  be 
obtained. 
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APPENDIX  A 


ADDITIONAL  OPTICAL  RESULTS 


For  reference  additional  optical  results  from  the  tank  experiments 
are  now  presented.  Spectral  albedos  and  extinction  coefficients  from 
the  growth  and  warming  phases  of  Experiments  3-6  are  shown  in 
Figures  A.l  -  A. 9.  In  each  figure  the  curves  are  numbered  by  site. 

Data  concerning  the  sites  displayed  in  the  figures  are  summarized  in 
Table  A.l. 

The  results  from  the  warming  phase  of  Experiment  4  are  somewhat 
unusual  (Figure  A. 5).  When  the  air  temperature  was  raised  from  -20°C 
to  0°C  the  albedo  increased  at  short  wavelengths  (less  than  600  nm)  and 
decreased  at  long  wavelengths.  This  was  due  to  the  surface  conditions 
for  Site  10.  The  surface,  having  undergone  draining  and  flooding,  was 
very  crumbly  with  surface  water  present.  The  albedo  was  enhanced  by 
the  crumbly  surface,  but  was  depressed  at  longer  wavelengths  by  the 


surface  water. 


Table  A.l 


Summary  of  Sites  Displayed  In  Appendix  A 


Experi¬ 

ment 

Site 

Thickness 

(cm) 

Time 

(hr) 

Air  Temp 
(°C) 

Surface 
Temp  (°C) 

Bulk  a 
Cloudy  Sky 

3 

2 

4.3 

9.0 

-10 

-5.4 

.12 

3 

7.6 

25.5 

-10 

-7.1 

.21 

4 

11.7 

-10 

-8.5 

.29 

6 

21.1 

178.0 

-10 

-9.9 

.37 

7 

27.2 

249.5 

-10 

-10.0 

.42 

8 

28.4 

336.5 

+1 

-1.3 

.33 

9* 

26.7 

339.5 

+5 

-1.3 

.31 

4 

3 

4.6 

18.5 

-20 

-12.6 

.10 

4 

7.6 

28.5 

-20 

-14.0 

.21 

5 

11.7 

42.5 

-20 

-16.4 

.31 

6 

17.8 

69.5 

-20 

-16.9 

.40 

8 

25.7 

140.5 

-16 

-16.0 

.53 

10^ 

27.9 

163.5 

0 

-1.0 

.53 

5 

6 

23.9 

71.5 

-30 

-24.2 

.63 

7 

29.0 

154.0 

0 

-1.5 

.42 

6 

2 

3.8 

5.7 

-31 

-13.0 

.13 

3 

7.6 

11.7 

-32 

-20.0 

.31 

5 

17.5 

33.0 

-36 

-26.5 

.52 

6 

23.5 

54.0 

-36 

-28.7 

.86 

7 

26.0 

74.0 

-37 

-31.9 

.90 

*  i I 

.8  cm  of  surface  water,  surface  very  crumbly  with  water  present 
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1  albedos  for  growth  phase  of  Experiment  3. 


Figure  A-l:  Spectra 
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Figure  A-2:  Spectral  albedos  for  warming 
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Figure  A-5:  Spectral  albedos  for  wanning  phase  of  Experiment  4. 
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Figure  A-7:  Spectral  albedos  for  wanning  phase  of  Experiment  5. 
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Figure  A-8:  Spectral  exti 
wanning  phase  of  Experi 
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-30,  and  -37°C)  and  from  water  of  different  salinities  (0,  16,  and  31  °/00)  was 
investigated.  The  experiments  were  conducted  using  a  tank  system  which  was 
designed  to  grow  ice  similar  to  that  found  in  nature,  allow  an  accurate  determi 
nation  of  the  state  and  structure  of  the  ice,  and  permit  in  situ  optical  . _ 
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^measurements  to  be  made.  Measurements  of  incident,  reflected,  and  transmitted 
irradiances  were  used  in  conjunction  with  a  modified  Dunkle  and  Bevans  photo¬ 
metric  model  to  determine  spectral  albedos  and  extinction  coefficients.  The 
thin  cold  ice  of  these  experiments  had  albedos  which  were  comparable  to  the 
values  for  the  thicker  warmer  ice  examined  by  previous  researchers;  however, 
extinction  coefficients  were  1.5  to  15  times  greater.  Qualitative  relation¬ 
ships  between  the  optical  properties  and  the  physical  state  of  the  ice  were 
observed.  As  the  ice  temperature  decreased  (and  the  brine  volume  increased) 
both  albedo  and  extinction  coefficient  increased;  when  the  ice  temperature 
dropped  below  the  eutectic  point  the  increase  was  drastic.  The  dependence 
of  albedo  and  extinction  coefficient  on  the  brine  content  of  the  ice  was 
found  to  be  complex  with  both  the  brine  volume  and  its  vertical  distribution 
being  significant.  Variations  in  ice  salinity  over  the  range  A  °/oo  to 
14  °/00  did  not  influence  the  optical  properties.  Increases  in  albedo  and 
extinction  coefficient  were  primarily  a  result  of  changes  in  the  ice  structure 
which  enhanced  scattering. 

^  A  four  stream  discrete  ordinates  photometric  model  was  developed  to 
treat  the  case  of  a  floating  ice  slab.  The  model  included  both  anisotropic 
scattering  and  refraction  at  the  boundaries.  The  effects  on  albedo  and 
transmittance  of  variations  in  such  model  parameters  as  the  single  scattering 
albedo  and  the  phase  function  were  investigated.  Using  one  and  two  layer 
models,  theoretical  albedos  and  transmittances  were  compared  to  experimental 
values.  The  four  stream  model  was  further  applied  to  investigate  such 
topics  as  the  depth  dependence  of  ir radiance,  cases  of  direct  beam  incident 
radiation,  the  spectral  dependence  of  albedo  and  transmittance,  and  the 
effects  of  surface  layers.  It  was  suggested  that  future  laboratory  research 
be  oriented  towards  determination  of  single  scattering  albedos  and  phase 
functions  for  an  extensive  range  of  ice  types. 


